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INTRODUCTION 

During  the  last  two  decades  there  has  been  a  growing  concern  by  ■ 
the  government  and  the  public  over  the  criteria  and  methods  by  which 
decisions  about  land  use  are  made.  An  increase  in  environmental  con- 
sciousness as  well  as  a  recognition  of  past  misuses  of  our  natural 
resources  led  Congress  to  pass  the  National  Environmental  Policy  Act  of 
1970.  It  is  the  intention  of  the  act  "to  encourage  productive  and 
enjoyable  harmony  between  man  and  the  environment;  to  provide  efforts 
which  will  prevent  or  eliminate  damage  to  the  environment  and  biosphere 
and  to  stimulate  the  health  and  welfare  of  man;  to  enrich  the  under- 
standing of  the  ecological  systems  and  natural  resources  important  to 
the  nation..." 

The  achievement  of  these  lofty  goals  presents  a  difficult  task. 
If  concern  about  environmental  quality  is  to  act  as  a  constraint  on 
urban  growth,  then  there  is  a  need  to  develop  and  apply  sensitive  and 
thorough  approaches  to  gauging  the  potential  impacts  of  urbanization. 
One  such  approach  is  the  use  of  computerized  information  storage  and 
retrieval  systems. 

The  application  of  computer  technology  has  become  an  increasingly 
important  facet  of  environmental  planning.  Several  studies  have  been 
completed  and  many  are  underway  which  demonstrate  the  application  of 
computerized  information  processing  and  data  manipulation  to  environ- 
mental planning.  The  intention  of  this  study  is  two-fold.  First,  to 
present  a  methodology  of  developing  a  computer  based,  geo-information 
data  bank;  and  second,  to  develop  a  package  program  to  display  that 
information  in  the  form  of  two-dimensional  maps. 
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A  geo-information  data  bank  may  be  defined  as  a  collection  of 
Information  describing  an  area  of  the  earth's  surface.  The  condition 
which  separates  geographical  information  systems  from  other  types  Is 
the  requirement  that  the  data  be  referenced  in  a  manner  that  will  allow 
retrieval,  analysis,  and  display  on  spatial  criteria.  The  system  for 
processing  this  information  is  an  nrganired  set  of  steps  or  procedures, 
which  accumulates,  manipulates,  and  disseminates  the  information  in  such 
a  way  that  it  can  be  used  as  a  basis  for  rational  decision  making.  The 
system  could  be  completely  manual  such  as  a  collection  of  maps  on  which 
information  is  recorded,  or  it  could  consist  of  manual  operations 
assisted  by  a  computer  with  its  associated  hardware  and  software  com- 
ponents which  may  be  utilized  to  perform  operations  on  the  data  base. 
Considering  the  myriad  of  interactive  variables  the  planner  must  con- 
sider, it  does  not  seem  that  a  completely  manual  operation  will  continue 
to  be  practical.  This  is  particularly  true  as  computer  usage  becomes 
more  and  more  available  to  local  and  regional  planning  agencies. 
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5 
To  demonstrate  the  application  of  a  geo-1nformat1on  data  system 
a  study  area  of  approximately  4  to  5  square  miles  has  been  selected. 
The  approximate  site  location  in  relation  to  the  state  is  shown  in 
Figure  1.  The  location  includes  the  City  of  Manhattan  and  its  sur- 
rounding area.  A  more  detailed  delineation  of  the  site  1s  shown  1n 
Figure  2. 

The  site  includes  parts  of  two  counties,  Riley  and  Pottowatomie. 
It  also  includes  two  major  rivers,  the  Kansas  and  the  Blue,  both  have 
rather  broad  river  valleys.  Manhattan  was  originally  settled  at  the 
confluence  of  the  two  rivers  and  most  growth  has  occurred  toward  the 
upland  area  to  the  west.  Recent  trends  in  land  development  are  toward 
the  east  into  the  river  valley.  With  the  construction  of  Tuttle  Creek 
Reservoir  on  the  Blue  River  flooding  is  no  longer  an  eminent  danger. 
However,  many  conflicting  land  uses  have  arisen.  Urban  expansion  in 
the  form  of  residential,  commercial,  and  industrial  land  uses  are 
competing  for  what  is  considered  prime  agricultural  land.  By  means  of 
a  geo- information  data  system  it  may  be  possible  to  1)  measure  in 
quantifiable  terms  the  impact  of  expanding  urban  development,  2)  deter- 
mine the  suitability  of  selected  sites  for  development,  and  3)  to  offer 
alternative  locations  for  orderly  growth  and  development. 

The  site  is  rectangular  in  shape  and  is  a  factor  that  was  perhaps 
Ill-considered.  A  project  scale  geo-information  data  bank,  such  as  this, 
would  be  better  defined  by  its  natural  watershed  boundaries.  A  county, 
regional,  or  state  data  bank  would  be  defined  by  its  political  boundaries 
which  may  or  may  not  lie  in  a  straight  line. 
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STUDY  AREA  LOCATION 


FIGURE  1     STATE  MAP  OF  KANSAS  ILLUSTRATING  THE  LOCATION  OF  THE  STUDY  AREA 


ill 


pip 

i±li:;:SS::::v:: 


FIGURE  2 


THE  STUDY  AREA 
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Data  Reference  System 

By  definition,  information  contained  in  a  geo-information  data 
bank  must  be  spatailly  related.  That  is,,  it  must  answer  two  questions: 
where  is  the  location;  and  what  are  the  geographic  data.  One  solution 
to  answering  these  questions  is  the  use  of  the  data  cell  system.  The 
cell  is  the  basic  unit  from  which  data  is  extracted.  The  Universal 
Transverse  Mercator  (UTM)  grid  is  one  system  used  to  locate  each  cell. 

The  UTM  system,  developed  by  the  U.S.  Army,  is  a  rectangular  grid 
system  designed  for  world  use  between  80°  south  and  80°  north  latitude.  A 
rectangular  grid  is  accomplished  by  mathematically  projecting  positions 
from  the  spheroidal  shaped  earth  onto  an  imaginary  cylinder  which  has  its 
axis  parallel  to  the  plane  of  the  equator  and  which  intersects  the  earth 
1n  such  a  way  that  zones  6°  in  an  east-west  direction  are  created.  There 
are  60  zones  each  6°  wide,  created  by  30  positions  of  the  cylinder.  A 
cylinder  is  a  "developable"  surface.  That  is,  it  can  be  rolled  flat  with 
no  distortion  so  that  a  uniform  grid  can  be  formed.  A  line  of  longitude 
serves  as  the  central  meridian  of  a  zone.  Multiples  of  6  longitude 
divide  each  zone.  Divisions  are  made  in  a  north-south  direction  by  lines 
of  latitude  in  multiples  of  8°.  Thus,  quadrilaterals  are  created  which 
are  6°  east-west  by  8°  north-south  and  each  of  which  is  given  a  unique 
Identification  called  the  grid  zone  designation.  Each  of  these  zones  is 
further  divided  into  100,000  meter  squares  each  of  which  has  a  unique 
2-1etter  designation.  All  recent  USGS  maps  bear  UTM  grid  ticks  every 
1,000  meters.  Connecting  these  ticks  will  result  in  a  grid  of  1-kilometer 
squares.  Because  the  basis  of  the  USGS  map  is  a  polyconic  projection  the 
UTM  grid  may  appear  slightly  askew. 
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The  numerical  value  given  to  the  grid  lines  is  based  on  convenience. 
The  arbitrary  value  for  each  central  meridian  within  a  zone  is  500,000. 
The  value  for  the  equator  is  0  for  the  northern  hemisphere  and  10,000,000 
in  the  southern  hemisphere.     Thus  a  grid  Tabled  4,351  .OOOrnN  is  4,351,000 
meters  north  of  the  equator  and  one  Tabled  714,000mE  is  east  of  the 
central  meridian  by  214,000  meters. 

The  UTH  system  has  provided  the  basis  of  the  inventory  grid  from 
which  cells  of  an  appropriate  size  may  be  constructed.     They  may  be  one 
square  mile,  one  square  kilometer  or  some  fraction  thereof,  one  quarter, 
one  sixteenth,  etc.     Cell  size  can  be  adjusted  to  fit  the  needs  of  the 
inventory.     The  selection  of  the  cell  size  should  be  made  with  careful 
consideration  of  the  data  requirements.     It  appears  to  this  researcher, 
that  very  small   cell  sizes  will   result  in  a  high  level  of  accuracy  but 
at  the  cost  of  spending  an  inordinate  amount  of  time  extracting  data;  a 
task  which  becomes  quite  tedious  over  a  long  period  of  time.     Consid- 
eration should  also  be  given  to  the  complexity  of  the  study  area.     An 
urban  area  or  a  site  of  highly  disected  topography  may  require  very 
small  cell  sizes  (1-acre)  while  a  large  area  of  uniform  topography  and 
large  open  spaces  may  be  adequately  studied  using  relatively  large  cell 
sizes  (1 -kilometer  squares). 

Using  the  1-kilometer  UTM  grid  as  the  base,  a  cell  size  of  250 
meters  square  was  selected  for  the  pilot  study.     This  results  1n  16  cells 
for  each  square  kilometer,  an  area  equivalent  to  6  hectares  or  approxi- 
mately 15  acres.     The  data  cells  are  identified  by  the  numerical  value 
given  to  the  UTM  grid  lines  with  some  modification  for  convenience. 
Because  the  maximum  length  of  the  study  area  is  only  44  kilometers  the 
first  three  digits  and  the  last  digit  of  the    UTM    grid  number  were 
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dropped  resulting  in  a  four  digit  fnasber.  Each  cell  is  then  numerically 
designated  by  the  value  of  the  grid  lines  below  and  left  of  the  cell.  . 
The  cell  is  then  identified  by  an  eight  digit  number.  The  first  four 
digits  represent  the  row  number  and  the  last  four  digits  the  column 
number. 

The  study  area  as  delineated  contains  1,760  cells  on  the  basis 
of  16  cells  per  square  kilometer.  The  reference  system  assumes  a 
rectangular  shape.  However,  irregular  shaped  areas  may  be  delineated 
by  assigning  values  of  zero  to  cells  outside  the  study  area. 
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FIGURE  4  CELL  IDENTIFICATION  SYSTEM 


Each  cell  has  a  specific  identification  number  based  upon  the  UTM  grid 
system.  The  first  four  digits  designate  the  row  number  and  the  last 
four  digits  designate  the  column  number.  The  cells  delineated  1n  the 
above  figure  are  0875  3850,  1000  4425,  1025  4200. 
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Variable  List  Development 

If  the  data  bank  is  to  provide  a  meaningful  tool  by  which  environ- 
mental planning  decisions  can  be  made  then  the  system  must  be  more 
accurate,  more  objective,  and  more  comprehensive  than  systems  used  1n 
the  past.  The  initialization  of  such  a  system  requires  that  careful 
attention  be  given  to  the  development  of  the  variable  list  which  forms 
the  basis  of  the  data  bank.  The  identification  of  resource  variables  is 
not  a  simple  task.  The  initial  listing  of  resource  variables  must  be 
broad  enough  to  insure  coverage  of  all  appropriate  resource  categories, 
it  must  be  specific  enough  to  provide  for  discrete  analysis,  and  it  must 
be  flexible  to  provide  for  additional  variables  as  more  descriptive 
information  becomes  available. 

For  the  purpose  of  this  study  a  resource  variable  will  be  defined 
as  a  primary  resource  element  or  data  type  which  is  descriptive  of  some 
natural  or  cultural  feature  of  the  environment.  An  attribute  is  an 
inherent  characteristic  or  quality  of  a  specific  resource  variable.  The 
following  section  examines  that  list  of  natural  and  cultural  variables 
selected  for  the  study  area.  The  data  processed  includes  physiographic, 
pedologic,  hydrologic  factors  as  well  as  man-made  features  and  influences. 

Natural  Features 

Sensible  environmental  planning  should  take  into  account  such 
factors  as  drainage  patterns,  slope,  elevation,  and  vegetative  cover. 
Encroachment  by  man  upon  these  ecologically  fragile  areas  may  result  1n 
significant  disruptions  to  these  natural  systems.  It  1s  important  that 
good  land  planning  take  into  account  the  location  of  these  natural 
resources  and  protect  them  from  damage  if  possible. 
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Data  relating  to  four  natural  features  were  Included  1n  the  pilot 
study.     They  we  I  e  as  follows:     1)  topographic  elevation,  2)  slope, 
3)  vegetative  cover,  and  4)  surface  water  patterns.     A  discussion  of  the 
importance  of  each  of  these  features  follows. 

Topographic  Elevation 

Topgrephic  elevation  is  the  measure  in  feet  of  the  earth's 
height  above  some  datum  plane,  usually  mean  sea  level.     Elevation  helps 
define  the  climatic  range  within  the  study  area  as  well   as  identifying 
flood  prone  areas.     Vegetative  associations  vary  with  elevation  in 
response  to  changed  micro-climate,  exposure  and  soil  conditions. 

In  order  to  achieve  greater  flexibility  than  programs  developed 
in  the  past  it  was  decided  not  to  predetermine  the  elevational   levels 
within  the  study  area.     Rather,  the  data  base  should  contain  specific 
elevi'.tional  information  for  each  cell.     The  computer  could  then  map  the 
topographic  relief  within  ten  levels  between  the  lowest  elevation  and 
the  highest  elevation.     The  centroid  elevation  of  each  cell  was  extracted 
and  stored  in  the  data  bank  for  this  purpose. 

Slope 

Slope  may  be  defined  as  the  degree  to  which  a  given  stretch  of 
land  is  inclined  away  from  the  horizontal.  The  formula  used  In 
calculating  the  slope  (as  a  percent)  may  be  expressed  as  the  ratio  of 
the  change  in  vertical  elevation  to  the  horizontal  distance  along  a 
straight  line,  times  100?«.  The  Importance  of  slope  in  land  planning 
cannot  be  overestimated.  Natural  systems  are  disrupted  by  development 
practices  which  misuse  land  of  various  slopes.  These  misuses  are  major 
factors  affecting  soil  erosion  and  surface  runoff.  Slope  may  be 
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regarded  as  one  of  the  more  critical  determinates  affecting  development 
in  urbanizing  regions. 

The  selection  of  slope  intervals  is  not  an  arbitrary  decision. 
Appropriate  criteria  for  the  selection  of  slope  intervals  should  relate 
to  the  geology  and  topography  of  the  region  under  study,  as  well  as 
the  use  to  which  the  slope  data  will  be  put  in  the  analysis  phase.     The 
Soil  Conservation  Service  and  the  Agricultural  Stabilization  and  Con- 
servation Service  generally  use  the  following  criteria  for  classifying 
slopes  according  to  agricultural  suitability. 

0-5%  slope      Low  Limitation 
5  -  105!  slope    Moderate  Limitation 
10  -  155.  slope    High  Limitation 
15%+  slope  Severe  Limitation 

Because  of  the  abrupt  changes  in  slope  associated  with  the 
topography  of  the  site  and  because  the  analysis  phase  is  intended  to 
include  uses  other  than  agriculture,  a  slightly  modified  classification 
system  for  this  study  was  developed.     The  following  attributes  of  the 
slope  variable  were  selected: 

Slope  attributes        0  -  2% 
3  -  6% 
7  -  12% 
13-20% 
21%+ 
It  should  be  noted  that  these  attributes  are  not  qualified  as  to 
the  degree  of  suitability  but  may  vary  according  to  the  analysis.     For 
Instance,  1n  selecting  sites  where  drainage  is  an  important  consideration, 
such  as  recreational  areas,  slopes  of  three  to  six  percent  would  be  more 
suitable  than  slopes  of  zero  to  two  percent. 
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Vegetative  Cover 

Vegetative  cover  may  be  defined  as  the  predominant  plant  life 
within  a  given  area.  The  type,  quality,  end  density  of  vegetative 
cover  is  an  Important  consideration  in  good  land  planning.  It  is 
necessary  to  maintain  good  natural  vegetation  in  some  areas  to  provide 
maximum  water  infiltration  and  prevent  soil  erosion  problems  from 
developing.  Woodlands  serve  as  greenbelts  that  aid  in  controlling  air 
and  noise  pollution.  They  also  provide  natural  screening  between  dif- 
ferent types  of  land  use  within  an  urban  environment.  Additionally, 
wooded  sites  provide  habitat  areas  for  wildlife  as  well  as  park  and 
recreation  facilities  for  an  urban  population. 

The  following  classifications  of  vegetative  cover  were  selected 
for  the  land  use  study. 

Marshes 

Prairie  Grass 

Shrubs 

Lowland  Forest 

Upland  Forest 

Surface  Water  Patterns 

Surface  water  patterns  include  streams,  rivers,  lakes,  ponds,  and 
all  man-made  water  bodies.  Knowledge  of  surface  water  patterns  is 
essential  to  environmental  planning.  Development  adjacent  to  drainage 
systems  may  drastically  alter  surface  runoff  which  affects  the  stream 
composition,  increases  water  pollution  and  has  an  exponential  effect  on 
higher  order  streams.  Additionally,  surface  water  provides  wildlife 
habitats  which  need  to  be  preserved.  Larger  lakes  and  rivers  provide 
sites  for  hiking  and  camping.  Surface  water  patterns  identified  in  the 
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pilot  study  include: 

Intermltent  Stream 

Creek 

River 

Lake  or  Pond  less  than  50  Acres 

Lake  or  Pond  greater  than  50  Acres 
The  classifications  of  Intermltent  stream,  creek,  and  river  are 
used  rather  than  a  system  of  stream  ordering  because  they  are  more 
readily  identified  on  USGS  maps.  The  distinction  between  lakes  and  ponds 
of  less  than  50  acres  is  the  basic  management  guideline  for  the  restric- 
tion In  the  use  of  power  boats.  It  should  be  mentioned  that  although 
aguifer  recharge  areas  were  not  included  in  this  study  they  represent 
a  critical  element  of  the  environment.  Development  on  these  ecolog- 
ically fragile  areas  may  drastically  alter  ground  water  and  stream 
elevations  as  well  as  affecting  the  quality  of  our  water  resources. 

Pedologic  Features 

Soils  information  plays  a  prominent  role  in  the  determination 
and  classification  of  land  for  various  purposes.  It  is  generally 
agreed  that  there  is  a  high  correlation  between  certain  soil  types 
and  agricultural  productivity.  Urbanization  has  added  the  need  to 
examine  soils  for  their  suitability  to  septic  tank  locations,  founda- 
tions for  buildings  as  well  as  roads,  sources  of  sand  and  gravel,  and 
other  developmental  phenomena. 

The  soil  classification  system  used  in  this  study  1s  that 
presented  In  the  Pottowatomie-Riley  Counties  Water  and  Sewer  Plan 
prepared  by  the  Oblinger  and  Smith  Corporation.  The  information 
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contained  1n  the  data  bank  pertaining  to  soils  has  been  extracted  from 

the  Soil  Structure  Map  contained  in  the  above  report.  The  following 
paragraphs  contain  a  detailed  discussion  of  the  various  land  use  classes 
and  their  capabilities. 

Class  1 

Lands  in  the  class  1  category  are  those  soils  which  are  best 
suited  to  cultivation  and  require  no  special  practices  for  erosion 
control  and  fertility  maintenance.  This  class  of  land  includes  soils 
which  are  deep,  friable,  and  silty  to  clayey.  These  soils  are  smooth  and 
usually  consist  of  well  drained  bottom  lands.  Also  included  in  this 
group  are  the  smooth  and  well  drained  stream  terraces  and  second  bottoms 
along  the  major  streams  and  valleys  of  the  study  area.  These  areas 
usually  include  the  best  cropland  and  should  not  be  compromised  to  the 
needs  of  expanding  urbanization. 

Class  2 

Soils  in  the  class  2  category  are  composed  mostly  of  sandy  soils. 
These  well  drained  soils  consist  of  alluvial  deposits  and  range  from  20 
to  38  inches  in  depth.  These  soils  are  moderately  fertile  and  suitable 
for  continuous  cultivation  requiring  only  moderate  agricultural  practices 
for  erosion  control  and  fertility  maintenance.  Wind  erosion  may  occur 
when  the  vegetative  cover  is  removed.  Of  the  5  solid  classes  identified, 
this  class  is  the  most  suitable  for  septic  tank  systems. 

Class  3 

Tight  clay  and  claypan  soils  compose  the  major  soil  group  in 
class  3.  These  soil  types  are  generally  found  on  slopes  of  three  to 
five  percent,  which  cause  the  surface  soils  to  have  less  depth 
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(8  to  10  Inches).  The  subsoils  are  composed  of  heavy  clays  to  silty  clay 
loams  which  extend  to  depths  of  20  Inches.  The  subsoil  is  sticky  or 
plastic  when  wet  and  very  hard  when  dry.  Class  3  soils  require  intensive 
or  complex  agricultural  practices  to  control  sheet  and  gully  erosion  and 
to  maintain  fertility.  These  areas  also  have  very  low  engineering  capa- 
bilities. 

Class  4 

Class  4  soils  consist  of  sandy  soils,  silty  to  clayey  soils  and 
tight  clay  and  claypan  soils.  The  topsoil  layer  is  thin  (8  to  10  inches) 
and  the  subsoil,  which  consists  of  silty  clay  to  clay  loam,  extends  to  a 
depth  of  three  feet.  These  soils  are  not  suitable  for  continuous  culti- 
vation but  may  provide  limited  production  of  cultivated  crops  with 
extensive  land  management  practices. 

Class  5 

Land  areas  in  class  5  are  typified  by  steep  slopes,  shallow, 
gravelly  or  stoney  soils,  severely  eroded  areas,  loose  sands,  and  fre- 
quently flooded  or  saline  bottom  lands.  These  soils  have  developed  over 
limestone  and  limey  shales  generally  on  steep  slopes  of  up  to  14  percent. 
The  topsoil  is  granular,  silt  loam  to  silty  clay  three  to  five  inches 
thick.  The  subsoils  seldom  go  deeper  than  sixteen  inches.  These  areas 
are  highly  susceptible  to  erosion.  They  are  unsuitable  for  cultivation 
but  are  suitable  for  grazing  or  woodland  and  limited  development. 

Cultural  Features 

The  geo-information  data  bank  should  also  contain  data  relating 

to  the  man-made  influences  upon  the  environment.  Three  cultural 

variables  were  selected  for  the  study.  They  were  as  follows: 
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1)  Activity  Patterns,  2)  Transportation  and  Utilities,  and  3)  Extrinsic 
Cultural  or  Historic  Features.  The  following  is  a  discussion  of  the 
importance  of  each  of  the  identified  features. 

Activity  Patterns 

The  activity  patterns  of  the  site  identifies  the  utilization  of 
the  land  by  man.  Potential  imbalances  created  by  human  activities  are 
reflected  in  present  land  use  patterns,  such  as  urban  development  in 
flood  plain  areas,  encroachment  on  prime  agricultural  lands,  and  the 
visual  pollution  along  transportation  corridors.  These  misuses  all 
impact  on  a  variety  of  natural  systems  including  wildlife,  the  soils 
system,  and  the  hydro! ogic  system.  Good  land  use  planning  should 
insure  a  proper  balance  between  man  and  nature. 

In  this  study  fourteen  attributes  relating  to  activity  patterns 
were  identified  and  extracted.  They  include: 
Residential 

Rural  (Agriculturally  related  housing) 
Suburban  (Low  density,  large  lot  sizes) 
Urban  (Medium  to  high  density,  fine  grain) 
Trailer  Homes  (Courts) 

Vacation  (Immediate  access  to  recreational  area) 
Commercial 
Industrial 
Agriculture 

Cropland 
Grazing  Land 


?/>■ 


Recreation 

State  and  Local  Park 

Golf  Course 

River  and  Lake  Zoning 
Institutional  Land 
Cemetaries 

Transportation 

Although  transportation  may  be  considered  as  a  land  use  activity, 
for  the  purpose  of  this  study  it  is  considered  a  separate  activity. 
The  transportation  network  is  one  of  the  primary  determinants  of  the 
direction  of  urban  expansion.  Residential,  commercial,  and  industrial 
land  uses  are  all  linked  directly  to  the  accessibility  of  transportation. 
The  location  of  existing  and  future  transportation  systems  is  related  to 
a  variety  of  natural  features  including  slope,  soil  type,  drainage 
patterns  and  vegetative  cover  as  well  as  land  use  activities. 

The  following  attributes  relating  to  transportation  systems  within 
the  site  were  identified. 

1.  Local  Serving  Street 

2.  Major  Arterial 

3.  County  Road 

4.  Unimproved  Road 

5.  State  Highway  (2-lane) 

6.  State  Highway  (4-lane,  divided) 

7.  Railway 


r 


LJa. 


3 


......i 


25 

Method  of  Extraction 

Having  determined  what  variables  are  to  be  included  in  the 
Inventory,  the  next  consideration  is  the  means  by  which  the  information 
will  be  obtained  and  stored  in  the  data  bank.  It  appears  that  there 
are  four  types  of  information  that  may  be  extracted.  Data  may  be 
extracted  by  name,  by  percent  coverage,  by  frequency  or  number  per 
cell,  or  by  presence  or  absence  within  each  cell. 

The  method  of  extracting  data  from  individual  cells  is  depen- 
dent on  the  nature  of  the  variable  and  the  use  to  which  the  information 
is  to  be  put.  Specific  natural  and  cultural  variables  were  extracted 
as  the  percentage  of  cell  that  a  particular  attribute  of  a  variable 
covers.  By  using  this  method  the  entire  coverage  of  a  cell  may  be 
stored  in  the  data  bank  rather  than  only  the  dominate  usage.  In  the 
case  where  more  than  one  attribute  is  presented  in  a  given  cell,  this 
system  becomes  much  more  discrete.  Other  features,  such  as  trans- 
portation Vines,  were  extracted  on  the  basis  of  number  per  cell. 

Finally,  certain  features  or  resources  such  as  hydrologic 
resources  were  extracted  by  name  (Z-digit  code)  on  the  basis  of  their 
presence  within  a  cell. 


METHOD    OF    EXTRACTION 


Source:  Miller,  1972 
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Hydro logic  System 
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The  Map  Programs 

The  Environmental  Resource  Planning  (ERP)  package  consists  of 
four  MAP  programs  written  in  PL1  ■  Each  1s  designed  to  produce 
spatially  related  maps  of  variables  that  have  been  extracted  using 
the  cell  system  described  previously.  Each  MAP  program  is  designed 
to  accept  one  of  four  types  of  data. 

MAPI  is  designed  to  accept  four  digit,  integer  data.  It  will 
produce  a  two-dimensional  map  in  ten  levels  ranging  from  the  lowest 
value  to  the  highest  value.  Uses  of  the  program  include  mapping 
topographic  elevations  which  may  include  recording  centroid  elevation, 
lowest  elevation,  and/or  highest  elevation  of  each  cell;  mapping 
density  characteristics  such  as  number  of  dwelling  units  per  cell; 
or  descriptive  characteristics  such  as  average  assessed  valuation  for 

each  cell. 

MAP2  is  designed  to  accept  data  stored  as  a  percentage  of  each 
cell.  It  can  map  up  to  fourteen  attributes  of  any  one  variable.  The 
program  produces  a  set  of  two-dimensional  maps  (one  for  each  attribute) 
illustrating  the  percentage  of  each  cell  occupied  by  the  attribute. 
Percentages  are  mapped  in  ten  levels  between  0  and  99.  Additionally 
it  gives  the  frequency  distribution  of  each  level  and  the  percentage 
of  usage  for  each  attribute  over  the  entire  site. 

MAP  3  produces  a  two-dimensional  map  of  up  to  ten  attributes 
of  a  single  variable  which  have  been  extracted  on  the  basis  of  predom- 
inate type.  That  is  the  attribute  which  has  the  largest  amount  of 
coverage  in  a  given  cell.  The  attributes  being  extracted  are  assigned 
a  coded  value  between  1  and  10.  Each  value  is  represented  by  an 
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overprint  character  on  the  computer  printout.  The  program  also  has  a 
provision  for  printing  out  the  appropriate  legend  and  the  frequency 
distribution  for  each  attribute. 

MAP4  is  designed  to  map  data  extracted  according  to  the  number 
of  occurrences  per  cell.  Unlike  MAPI  this  program  does  not  require  a 
data  card  for  each  cell,  however,  it  will  only  accept  integer  values 
between  1  and  10.  As  many  as  35  attributes  of  a  variable  may  be 
stored  in  a  single  data  file.  The  program  will  produce  a  set  of  two- 
dimensional  maps  (one  for  each  attribute)  illustrating  the  density 
pattern  of  each  resource  and  the  frequency  distribution  for  each  value. 

In  certain  cases  data  that  is  extracted  for  use  in  one  program 
can  be  converted  for  use  in  another  program.  For  instance,  Appendix  B 
includes  a  program  that  will  analyze  a  data  deck  in  which  a  variable 
has  been  extracted  by  percent  of  cell  and  produce  another  data  deck 
in  which  the  information  is  converted  to  predominate  type. 
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MAPI:  PROCEDURE  CPTICNS  (MAIN); 

/a*************** *************#<•****>«*****************  *************<**/ 

/#       MAPI  IS  A  PROGRAM  OESIGNED  TO  MAP  SPATIALLY  RELATED  INTERVAL 

*/ 

/*  DATA.    IT  TS  MOST  USEFUL  IN  MAPPING  DENSITY  PATTERNS  OR  CONTOUR 

*/ 

/*  INTERVALS.   THE  DATA  WILL  RANGE  IN  VALUE  FRO"  A  MINIMUM  TO  A 

*/ 

/*  MAXIMUM  IN  UNIFORM  INCREMENTS.   THE  MINIMUM  THAT  MAY  BE  USED  IS 

*/ 

/*  ZERO  AND  TFE  MAXIMUM  IS  9999.99.   THE.  USER  MUST  PROVIDE  A  NUMBER 

*/ 

/*  OF  CARDS   IN  THE  SOURCE  DECK.   THEY  ARE: 

*/ 

/*            IMAX  -     (NUMBER  OF  ROWS)  ; 

*/ 

/*            JfAX  -     (NUMBER  CF  COLUMNS)  ! 

*/ 

/*            KMAX  =  (NUMBER  CF  VARIABLES  TO  BE  MAPPED)  ; 

*/ 

/*            CFIL-NUMB  =  (TOTAL  NUMBER  OF  CELL  IN  STUDY  AREA); 

*/ 

/*            MINI*)  «  (LOWEST  VALUF  OF  EACH  VARIABLE  TO  BF  MAPPED)  ; 

*/ 

/*            MAX!*)  =     (HIGHEST  VALUE  OF  EACH  VARIABLE  TO  BE  MAI'PECl? 

*/ 
*/ 

/* 

/*       ADDITIONAL  CARDS  FOR  PRINTING  OUT  APPROPRIATE  LABELS  MUST 

*/ 

/*  ALSO  BE  SUPPLIED  BY  THE  USER.   THE  VARIABLE  NAME  'TITLE'  IS  A 

*/ 

/*  CHARACTER  STRING  OF  50  CHARACTERS  USED  TO  IDENTIFY  THF  VARIABLE 

*/ 

/*  NAME.   'LABEL!*)'  IS  A  VARIABLE  NAME   DESCRIBING  THE  SPECIFIC 

-■'.' 

/*  ATTRIBUTE  OF  THE  VARIABLE  BEING  MAPPED.   THERF  IS  ONE  LABEL  FOP 

*/ 

/*  EACH  ATTRIBUTE.   THE  MAXIMUM  LENGTH  IS  100  CHARACTERS. 

*/ 

/*  THE  USER  CARD  SHOULD  FOLLOW  THE  FORMAT  BELOW: 

*> 

/*            TITLE  =  'NAME  OF  VARIABLE'; 

*J 

/*            LABEL!*)  =  'NAME  OF  ATTRIBUTE'; 

*/ 

/*  THE  ASTERISK  (*l  REPRESENTS  TFE  SEQUENTIAL  NUMBER  OF  EACH 

*/ 

/*  ATTRIBUTE. 

*/ 

/*      THF  USER  MAY  ALSO  USE  UP  TO  FIVE  LINES  ON  THE  COMPUTER 

*/ 

/*  PRINTOUT  TO  PROVIDE  INFORMATION  PERTINENT  TO  HIS  SPECIAL  PROJECT 

*/ 

/*  OR  STUDY.   THE  FIRST  LINE  MAY  CONTAIN  UP  TO  100  CHARACTERS  AND 

*, 

/*  THE  REMAINING  FOUR  I.INFS  CCNTAIN  50  CHARACTERS  EACH.   IF  THF  USER 

*i 

/*  DOES  NOT  WISH  TO  USE  THE  LINE  OPTION  HE  MUST  USE  THE  STATEMENT: 

*l 

/*           LINEU  )  ^  ■  ' ; 

*i 

/* 

*t 

/#***# ********************************************************* *******  i 

DCL  (MIN(I),  MAX(l))  FIXED  DEC  (7,2); 

CCL  CELL_NUMfl  FIXED  DEC ; 

CCL  LINE15)  CHAR  (50); 

DCL  TITLE  CHAR  (50) 5 

CCL  LABEL!  1  '< )  CHAR  (  100)  ; 

/*************«4*************>;  ****************************************  J 

/*  USER  CARDS  FOLLCW  */ 

IMAX  = 

JMAX  = 

KMAX  = 

CELL_NUMB  = 

MINI*)  =  USER  SUPPLIED  INFORMATION 

MAX!*)  = 

TITI.F  = 

LABEL!*)  = 

LINE  I* I  = 

/********************************************************************* 

begin; 
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OPEN  MIT  (SYSPRINT)  PAOESI 1 1  (  00  ) : 

DCL  MAPI [MAX, JMAX.KMAa)  FIXED  BIN; 

DCL  NUMPI IMAX.JMAX.KMUX)   FIXED  DEC  (7,2); 

DCL  (ROWtlMAXI,  COL(J"AX)l  FIXEO  HIM; 

DCL  LEVEL(KMAX.IO)  FIXED  DEC; 

DCL  (OIF,  INCRI  FIXED  CEC; 

DCL  FREC(KWAX,10)  FIXED  DEC  INITIAL  ( ( KMAX* 10 ) 0  )  ; 

OCL  PERCENTtlOl  FIXED  DEC  (7,21; 


/*  DETERMINE  THE  DIFFERENCE  ('DIEM  BETWEEN  THF  MIN(*I  AND  MAX(*I     */ 
/*  VALUES  AND  DIVIDE  BY  1C.   THE  RESULTANT  IS  THE  INCREMENT  C1NCR')  */ 


/*  BY  WHICH  TFE  DATA  WILL  BE  MAPPED  IN  10  LEVELS. 


*/ 


DO  K  =  1  TO  KMAX; 

DIF  =  MAX(K)  -  MINIK); 

I  NCR  =  DIF/IO; 

LEVEL!  K,  I  1  =  HMKI  ♦  INCP; 

DC  N  =  2  TC  10; 

LEVEKK.N)  =  LEVELIK.N-1)  ♦  INCR; 
ENC; 
END; 


/*  READ  IN  VALUES  OF  THE  EXTRACTED  DATA  SO  A  MAP  VALUE  MAY  BE 
/*  ASSIGNED  FOR  EACH  CELL. 


*/ 

*> 


DC  1  -  I  TC  IMAX; 

DO  J  =  I  TC  JMAX; 

GET  EDIT  (ROW(I),  COL(J))  (CCL(l),  F  1  <. )  ,  C0LI6J,  F(4)l; 
DC  K  =  1  TC  KMAX; 

GET  EDIT  (NUMBU  ,J,K)  )  (F(4,0)  );    /***CHANGE***/ 
CO  K  =  1  TC  io; 

IF  NUfB(I,J,KI  <=  LEVEL(K.N)  THEN  00; 
MAPd  ,  J,K)  »  N; 
FRECtK.M  =  FREQ(K,N)  +  1; 
GC  TO  MATCH; 

END; 

END; 

MATCH:       ENO;  /*   K  LUOP   */ 

END;       /*  J  LOOP   */ 

END;        /*   I  LOCP   */ 

/*  AFTER  ALL  THE  DATA  HAS  BEEN  TRANSFORMED  TO  MAP  VALUES  OVERPRINT    */ 
/*  CHARACTERS  ARE  ASSIGNFD  TO  EACH  VALUE  TO  DISPLAY  THE  MAP  PRINTOUT  *> 

/a********************************************************************/ 


DCL  SYM(IO)  CHAR11I 
SYMt  1  !-•  .• 
SYM!2)=' I' 
SYMt  ))='-' 
SYP(4)='  ♦■ 
SYM(5>=' X' : 
SYMI  6)=«C 
SYM(7)='=' 
SYMt  G)  =  <  | ' 
SYM(9)»«X«  ; 
SYMt  1C)  =  M'  ; 
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DCL    ARRAYCt.JMAXI    FIXED    BIN    INITIAL     !K*JMAXSO); 
DCL    I  ITEVP,NFIX(<il  )    FIXED    WIN  ; 

DO   N    =    I    TO   JNAX    BY    2i 
ITEMP    =    CDL  IN)  ; 

on  i  =  i  to  3; 

NFIXl  I  )    =    I  TEMP/lO**K-I  I  ; 

ITEMP  =  ITEMP-INFIX( I I*10**(4-I )); 

ARRAYU  ,NI  =  NFIXl  I  )  ; 

END! 
ARRAYS,  N)  =  ITEMP; 
E  ND ; 

/**#*** ***<***  ***********  « *********** ***********************  **********/ 


DO  K  =  I  TO  KMAX; 

PUT  EDIT  (T  1TLE)  ( A(5C)  ) ; 
DO  K  =  I  TG  4; 

PUT  SKIP  EDIT  ['  ■)  (C0L16),  AID); 
DC  N  =  1  TC  JMAX  RY  2; 

PLT  EDIT  (ARRAY(M,N>)  (F(2)l; 
END; 
END; 
PUT  EDIT  (!'*'  DO  I  =  b    TO  JMAX+9))  (CCL(I)t  A(ll); 
PUT  ECIT  (•*',  '*')  (C0LI6),  All),  COL  (  J  MA  X  +  <3 )  ,  All)); 
DC  I  =  1  TO  IMAX; 

PUT  EDIT  IR01-.II),  '*')  (CCLll),  F[4),  X(l),  All)); 
DO  J  =  1  TO  JMAX; 

IF  MAPI  I  , J,K)  >  5  THEN  DOS 

PUT  EDIT  ISYM(6>)  lCOLU+7),  4(1)1; 

end; 
no  N  -  1  TO  5; 

IF  MAPI  I , J, K)  =  N  THEN  DO; 

PUT  EDIT  (SYK(NI)  (C0LU*7I,  All)); 
END; 
END; 
END; 
PUT  SKIPIO); 
DC  N  =  7  TC  10; 

DC  J  =  I  TO  JMAX: 

IT  MAPI  I, J, K)  >=  N  THEN  DO; 

PLT  EDIT  (SYM(NI)  IC0LU+7),  All)); 

end; 
end; 
put  skipio) ; 
end; 
put  edit   [•   *   ',   row(i))    (coll  jmax*8>  ,  a(3),  f  ('.)); 

END; 
PUT  EDIT  (•*',  **•)  IC0L16),  All),  C0LIJMAX+9),  All)); 
PUT  EDIT  ((•*'  DO  I  =  0  TO  JMAX*9))  (COLII),  All)); 

DC  »  =  1  TC  '. ; 

PUT  SKIP  ECIT  (•  ')  (C0L16),  & (  1 1 )  ! 
OC  N  =  1  TO  JMAX  RY  2; 

PLT  EDIT  IARRAY(M,N)I  (  F  (  2  )  )  ; 

END; 
END; 
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/*  THE  REMAINING  PORTION  OF  TFF  PROGRAM  PRINTS  THE  LFGFND  AND 
/*  OTHER  INFORMATION  PROVIDED  BY  THE  USER  BELOW  THE  MAP. 


*/ 


PUT  SKIP!?)  EDIT  (IA3FLIK))  lAttOOt  t  i 
IF  L INF  I  1 )-=  '  '  TFPN  DO; 
PUT  SKIP!?  I  ; 

PUT  EDIT  I  LINE! I)  I  (A(IOO) ) : 
DO  N  >  2  TO  5i 

PUT  EDIT  (LINF(N)l  (COL(50),  A(50)>: 
END; 
END; 

PUT  SKIPI3)  EDIT  CCATA  MAPPFD  IN  10  LEVELS  FROM  THE  MINIMUM  VAI.L 
F*,  MIN(K),  'TO  THE  HIGHEST  VALUE',  MAXIK))  IA147),  F(7,2>,  Xll), 
A (22) ,  F(7,2) I; 

PUT  SKIP13)  EDIT  ('VALUE  RANGE  APPLYING  TD  EACH  CELL')  (A(60)l: 
PUT  EDIT  ('MINIMUM*,  MIN1K))  IC0LI2),  A(10),  C0LI2  1I,  F(7,2)); 
DC  K    =  1  TC  9! 

PUT  EDIT  (LEVEL(K,NH-.OL)  (X(3),  F(7,2)>; 

end; 
put  edit  ('maximum')  (c0l(2),  a1161); 

DO  N  =  1  TO  10; 

PUT  EDIT  (LEVFL(K,N>)  (X(3),  F(7,2)l; 
END; 
PUT  SKIP(2)  EDIT  ('PERCENTAGE  OF  STUDY  AREA  INCLUDED  AT  EACH  LCVEL'I 

(AUCOI  >  ; 
PUT  SKIP  EDIT  ('  ')  (CCLI17),  A(l)); 
DO  N  =  I  TO  10; 

PERCENT(N)  =  ( FPFQt K, N ) /CELL_NUMB)* 100 : 
PUT  EDIT  (PERCENT(N)I  (  X  (  4  )  ,  F(6,2I); 

END; 
put  skip12)  edit  cfrecufncy  distribution  of  data  cell  values  in  ea( 
h  level' )  !a( 100) i ; 

put  skip12)  edit  ('levels      0')  (x(l),  aihii; 
do  nlfvel  =  1  to  10; 
put  edit  (nlevel)  ( x ( 9  i  ,  f(2) )  ; 
end; 
put  skipu  )  ; 

PUT    EDIT     (('='     DC     I     =     10    TO    120))     (COl(I),     A(ll); 
DO    X    =    1    TO    5; 
IF    X    =    3    THEN 
DC; 

PUT    SKIP    EDIT     ('SYMBOLS')     (X(l),     A17M; 
END; 
ELSE    DC; 

PUT  skip; 
END; 

PUT  EDIT  (■ I  II  Ml!  I! +4-4-H-  +  +  +  +  XXXXXXXXX  OOOOi 

0000  COOOOOOOO  OOOOOCOOO  COOOOOOOO  000000000')  (C0L121),  A(IOO)); 
PUT  SKIP  (0); 
PUT  EDIT  («=««—«  ===*.==«=  =  =  =  =  =  =  =  =  =  »»»==»===t)  (C0H81I, 

A  ( 4  C 1  )  i 

PUT  SKIP  !0) ; 

PUT    EDIT     enillllll     lllllllll      lllllllll'l     (C0L<91),    A(30)l; 

put  skip  lo); 

PUT  EDIT  ('XXXXXXXXX  XXXXXXXXX')  (CCL  (101),  A(20)); 
PUT  SKIP  (0) ; 
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PU1  EDIT  ('**## «««»#•>  (COL (111)  i  A(10)); 

END; 

PUT  SKIP; 

PUT  EDIT  ((•  =  '  DO  I  =  II)  TO  1HO))  (COL(I),  AilDi 
PUT  FOIT  (•FREOUeNCV'l  (CGL12),  AU61); 
DO  N  =  1  TO  10; 

PUT  rr.IT  (FREOIK.N))  (X('«),  F(6)l; 
END: 
END;         t*       K  LCOP   */ 
FND!      /*  BEGIN  8LCCK   */ 
END  MAPI; 
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MAP2:  PROCEDURE  OPTIONS  l«MN)l 

/*  SITEMAP  II 

/* 

/* 


BY  THOMAS  M.  KUNTZ 

DIPT.  OF  REG.  f.  COMM,  PLANNING 
KANSAS  state;  UNIVERSITY 


/*> 

/* 

/* 

/* 

/* 

/* 

/* 


*******  ** 


,*** ***  **  *************** *** ** ** *****<  *******  *******  < 


MAP?  IS  A  SPECIALIZED  PROGRAM  'IF  THE  SITEMAP  II  PACKAGE  AND 
IS  DESIGNED  TO  MAP  SPATIALLY  RELATED  DAT/  THAT  IS  FXT'ACT^D  AS  A 
PERCENT  TF  FACH  CFLI  .   THE  DATA  MAY  RANGE  IN  VALUE  FRO"  ZFRO  TO 
1  00.   AS  MANY  AS  K  ATTRIBUTES  [IF  FACH  VARIABLE  MAY  HE  MAPPED. 
THE  FUILCWIKG  CARDS  ARF  REQUIRED  F*Y  THE  USER  T0  ESTABLISH  THE 


/*  BOUNDARIES  OF  THE  SITE  OR  STUDY  APFA. 


/* 

/* 
/* 
/* 
/* 
/* 
/* 
/* 
/* 
/* 
/* 
/* 
/* 
/* 
/* 
/* 


I«AX  "  (NIMBFR  OF  ROWS) ; 
J«AX  =  (NUMBER  OF  COLUMNS  I  I 
KMAX  ■-  (NUMBER  OF  ATTRIbUTES 
CELL  NUMB  =  (TOTAL  NUMBER  OF 


FUR  VARIABLE) 
CELL  IN  STUDY 


AREA); 
ITLES  AMD 


ADDITIONAL  CARDS  c CR  PRINTING  OUT  APPROPRIATE 
LEGEND  MUST  ALSO  BR  SUPPLIED  BY  THE  USER.   THC  VARIABLL  NAME 
•TITLE'  IS  A  CHARACTER  STRING  OF  50  CHARACTERS  USED  TO  IDENTIFY 
THE  VARIABLE  NAME.   'LABEL!*)'  IS  A  VARIABLE  NAME  DESCRIBING  THE 
SPFCIFIC  ATTRIBUTE  BEING  MAPPED.   THERE  IS  ONE  LABEL  CAF  D  FOR 
EACH  ATTRIBUTE  (MAXIMUN  IS  141  AMD  THE  LABEL  IS  LIMITED  TO  100 
CHARACTERS. 

TITLE  =  'NAME  OF  VARIAEl q ' ! 

L  ABEL ( *  >  =  'NAME  OF  ATTRIBUTE'; 
THE   ASTERISK  REPRESENTS  THE  SEQUENTIAL  NUMBER  OF  FACH  VARIABLE 

AT1F I  BUTE. 

THE  USER  MAY  ALSO  USE  LP  TO  FIVE  LINES  OK  THE  COMPUTER 

TO  PROVIDE  INFORMATION  PERTINENT  Ti]  HIS  SPECIAL  PROJECT 
OR  STUDY    THE  FIRST  LINE  MAY  CCNTAIN  UP  TO  11(1  CHARACTERS  AND 

GAINING  FOUR  LINES  CCNTAIN  50  CHARACTERS  FACH.   IF  THF  USFR 
T  WISH  TO  LSF  TF-F  LINE  OPTION  HE  MUST  USE  THF  STATEMENT! 


THF.  RE> 
DOES  NOT 


/* 

/*  PRINTOUT 

/*• 

/< 

/  + 

f*  LINF(l)  =  '  '■ 

/* 

/***«  ****************** **#< 

DCLUMAX,  JMAX,  KMAX)  F-IXFl. 
DEL  C  El.L_NUNB  F  IXFD  BIN; 
DCL  T  ITLE  CHAR  1  VI)  ; 
TCL  LABEL ( 14)  CHAR  ( 100) : 
DCL  LINE  (5)  CHAR  (100)  ; 

/ ****rt *********  ************ 

/*    USER    CAFES    FCl  LOW    */ 


BIN; 


|c*  ****#***************#*■***#  ******** 


*/ 
*/ 
*/ 

<*/ 
*/ 
*/ 
*/ 
*/ 

V 

*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
■>/ 
*/ 

**/ 


>**** *** ***************** ***#***#**/ 


I  MAX  = 
JMAX  = 
KMAX  = 
CELL.NUMB 

TITLE  = 

la  ha  (*)  = 

LI  NEC:- I  = 


/*  USER  SUPPLIED  INFORMATION  */ 


/****. 


*********•>**********»********>!  *************#***<  ***************> 
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BEGIN; 

OPEN    FILE     (SYSPRINTI    PAGFSIZE!  160»  5 

OCL    ATTR1KM.AX)     FIXED    Ht"C(5,2Ji 

DCL     (RO«(I"AXI,     C?HJ"AXII     FJXEf)    BIN; 

DCL    MAPI  I«AX,.IMAX,KMAX)    FIXED    DECC5,2); 

OCL    REMAP! IMAX, JMAX.KMAX)     FIXED    BIN; 

PCL    FREQIKMAXtll)     FIXED    PIN    I  N  I  T  I  A.L     (.  (  KMAX*  11  I  0  )  : 

DCL    PERCENT (KMAX)    FIXFD    DEC    (10,2)     INITIAL    l(K«AX)OI 

DCL    JUMP    CHAR    (  1  I  ; 


INPUT    DATA    AND    ASSIGN    VUUE    TO    MAP  (I,  J, K) 


*/ 


DO  1  =  I  Tf  [MAX: 
DO  J  =  1  TO  JMAX; 

GFT  EDIT  (RCWII).  COL(J),  JUMP)  (COLIll,  F(4),  X(ll 

a  ( 1 1 ) ; 

DO  K  =  I  TO  K m A  X ; 

GFT  FDIT  (ATTR <K1  )  (F<5,2))J 

"API  I  , J,K)  =  ATTKIK) ; 

PERCENT  (  K)= PERCENT  (K  H- ATTR  IK)  ; 


:  (  4  )  , 


/*    REASSIGN   CORRESPONDING    VALUE    TO    R E MAP  I  I , J , K I 


*/ 


SEEK: 


FIND: 


ADDER    «    0; 

kgunt  =   n; 

IF  MAPI  I, J,K]  = 
REMAP! I , J, K) 
FREO(K,KOUNT 

go   ro   find: 

END; 
ADDER  =  ADDER  + 
KCUNT  =  KCUNT  t 


/*  MAP  VALUE  COUNTER  '</ 
/■"    RFMAP  VALUE  COUNTER  */ 
ADDEP  THFN  00 ; 
=  KUUNT; 
-I)  =  FREOlK.KOUNT+D  +  l; 


lo; 


do; 


IF  MAPI  I. J, K)  <=  ADDER  THEN 
REMAP! I , J, K )  =  KOUNT ; 
FRFQtK,KOUNT+l)  =  FRFOIK, KCUNT+1 1+1 i 

go  to,  find; 
FND; 

IF  ADDER  ->=    1. 00  THFN  GO  TO  SEEK! 
END; 


END; 
END; 


/*»*##****  *#*i  ************** ***********#******#*******#* *******  -:&#****/ 
f*  AFTER    ALL    THE    EXTRACTED    DATA    VALUES    HAVE    flf'FN    CONVERTED    TO  */ 

/*    REMAP    VALUES    A    SPATIALLY    RELATED    MAO    USING    CHARACTER     SYMBOL <     IS         */ 
/*    PRINTFD    OUT.       THE    FOLLOWING    CHARACTER     SYMBOLS    A?H    USED!  */ 

/**  ******************  <>*********************************#******<  *******/ 


DCL  SY" 
SYMI  1  )  = 
SYP(2)= 
SYMI  3)  = 
SYMCt  1  = 
SYM(5)= 
SYMU>I  = 
SYM(7)= 

SYf i at  = 


110)     CHAR  (  I): 
I 
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SY«( 9) =  •  y : 

SYM(10>='*' i 


/*  ARRAY(*i*)  IS  A  TEMPORARY  STORAGE  SPACE  TO  ALLOW  THE  UTH  GRIE      */ 
/*  NUMBERS  USEC  FMR  COLUMN  REFERENCE  POINTS  TO  RE  PRINTED  VERTICALLY  *l 


DCL  APRAV  (4, JMAX)  FIXED  BIN  INITIAL  I  ('.*  J  MAX )  0  )  ; 
DCL  (  ITFMP,  NFIXUI  )  FIXCU  BIN? 


DC  N  =  1  TO  jr*AX  RV  ? \ 
ITFMP  =  COL  IN  I  : 
DO  I  =  I  TO  3 i 

NFIX(I)  ■  ITFMP/10**(4-1 ) ; 

ITEMP    =     ITEMP-INFIX!  I  )*10**(4-I  I  )  ; 

APRAYt! ,N1     =    NFIXI I) ! 

END; 
ARRAYC-..M     =     I  TEMP; 
END; 


/*  BEGIN  PRINTING  OUT  MAP  PICTURE 


*/ 


DO   K    =    no   K»AX; 

PUT    EDIT     (TITLE  I     ( A ! 50 )  )  ; 

PUT  skipu); 

DO  M  =  1  TO  4: 

PUT  SKIP  EDIT  (•  ')  (CCLI6I1  All))! 
DC  N  =  I'O  JMAX  RY  2; 

PUT  LDIT  (ARRAY(M,N)I  (F(2>); 
END; 
END; 
PUT  EDIT  (!•*•  DO  I  =  6  TO  JMAX  +  <* ))  !CCLd),  AUI): 
PUT  ECIT  {'*•.  '*')  (CCL(6),  A(l),  COL  ( J  MAX  +  <5  )  ,  Allll: 
DO   1=1  TO  IMAX; 

PUT  EDIT  (RCW(I)t  '*')  (COL(l)i  F  (4  I  ,  X(1),A(1I); 
DC  J  =  1  TO  JMAX; 

IE  REMAPd  ,J,K)  >  5  THEN  DO; 

PUT  EDIT  (SYMI6I)  (C0L(J+7)i  A(l)l; 

end; 

FLSE  DO; 

DC  N  =  1  TO  5; 

IF  REMAPd  ,  J,K)  =  N  THEN  DO; 

PUT  FDIT  (SYMINI)  (COLtJ+7),  Allll; 
END; 
END; 
END; 
END; 
PUT  SKtPtOI; 
DC  N  =  7  TO  10: 

DC  J  =  1  TO  JNAX; 

if  remapd  ,j,k>  >=  n  then  do; 

put  edit  (sym(mi)  (colu+71.  au)>; 
end; 
end: 
plt  skipioi; 
enc; 

PUT  EDIT  ('  *  'i  ROWdl)  «Cni.(JMAX  +  fi)f  A(3I.  E  ( ',  I  ) 
END; 
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PUT     EDIT     (•*",     •*')     !COL((>),     ftllli     COLIJMAX+Q),     A  (  I )  }  ; 
PUT     FDIT     ((•*■     DO     1     --     6    TP    JMAX+9))     (COM!  It    All)); 

DC   m   =    i    m   4; 

PUT    SKIP    EDIT     ('     •)     IC0LI6I,     Mllli 
DC    N    =    I    TIj    JKAX    BY    ?; 

PUT    LD I T     (ARRAY(M,N) )     ( F ( 2  I  I ; 

END; 

END; 

/*  THE  REMAINING  PORTION  OF  THE  PROGRAM  UTILIZES  THE  USER  */ 

I*    SUFPLIFD  INFORMATION  TC  LABEL  THE  MAP  AND  PROVIDE  THE  APPROPRIATE  */ 
/*  LEGEND.  */ 

PUT  SKIPI2)  EDIT  (LABFL(K))  (AllOOIIi 
PUT  SKIP12); 

IF  LINE!  1)  -=  ■  ■  THEN  DO; 

PUT    EDIT    ILINE(1 ) )     IAI10OII; 
DO    N    =    2    TO    5; 

PUT    EDIT     (LINF(N)I     (COL150),     »(50)l! 
FND; 
END; 
PUT     5KIPI3)     nDIT     ('DATA    IS     EXTRACTED    AS    A    PERCENT    OF    EACH    CELL    AN 
D'l    (At  100)); 

PUT    SKIP    EDIT     ('MAPPED     IN     10    LEVELS    BETWEEN    EXTREME    VALUES    OF     0.0 
0    AND     ICO. 00    PERCENT')     (A(100))i 

PUT    SKIPI3)     EDIT    ('VALUE    RANGE    APPLYING    TC    EACH    LEVEL')     (AtlOOl): 
DCL     (KIN,    MAX)     FIXED    DEC     (6,?); 
MIN    =    0.C1 ; 
MAX     =     100.00: 
ADDER    =    »U; 

PUT    EDIT    ('MINIMUM',    MIN)     (C0K2J,    A(IO),    COLI22),    F(6,2|); 
ADDl:    ADDER    =    ADDER    +    10. CO; 
IF    ADDER    <    MAX    THEN    DO; 

PUT    EDIT     (ADDER)     (  X  (  4 ) ,     F ( 6 , 2 ) I  ; 

GC    TC    ADOli 

END; 

PUT    ECIT     ('MAXIMUM')     <CCL(2),     A(16)l; 
ADDER    =    0; 
ADD2:     ADDER    =    ADDER     +    10.00; 

PUT    EDIT     (ADDER)     (XK),     F(6,21); 
IF    ADDER    <    MAX    THEN    GO    TO    ADD-?; 

/*    CALCULATE    PERCENTAGE    OF    TOTAL    SITE     IN    WHICH    RESOURCE     IS    PRFSFNT  */ 

DCL    PERCEI\T_T(H     FIXED    DEC     (6,?); 
PFRCENT_TOT    =     (  PERCENT  (K  )  /CEL  L_NUMP;  I  *100  ; 

PUT    SKIP(2)    EDIT     ('TOTAL    PERCENTAGE    OF    AREA    IN    WHICH    RESOURCE     IS     PRE 
SENT    FGUALS',     PERCENT_TOT  )     (A(6l),     F(6,2))S 


/*    PRINT    OUT    FREQUENCY    DISTRIBUTION 


*/ 


PUT    SKIIM2)     EDIT     ('FREQUENCY    DISTRIBUTION    OF    DATA    CFLL    VALUES     I.N    EAC 
H    LEV  El  •  )     I  A!  ICO)  )  ; 

PUT    SK1IU2I     EDIT     ('LEVELS  0')     (X(l),     A  (  1  't  1  )  ; 

00    NLEVEL    =     1    TO     10; 

PUT    LO  IT     (NLEVEL)     (X  (B  ),    r(?)  )  ; 


• 
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ENOi 
P 11  SKI P(l I ; 

PUT  EDIT  ((*  =  '  Of)  I  =  10  TO  120)1  (Ctll.ll).  41111! 

no  x  =  i  to  ^: 

IF  X  =  3  THFN 

CO; 

PUT  SKIP  EPIT  (•SYfBCLS')  (X(ll.  A(7)l; 

END; 
ELSE  CO; 

PUT  SKIP: 

FNn; 
puT  Enu  (1 Illllllll  +  +  +  +  *+  +++  xxxxxxxxx  ooooo 

OCOO  COCCCCOGO  CCOCCCCOO  CGCOOOOOO  0000000)0'!  (C0L121),  A ( 1 00  )  )  ; 
PUT  SKIP  (0); 
PUT  EDIT  [t=========  =========  =========  =========•)  (COL(Bl) , 

AI4CI ) ; 
PUT  SK!"  (0) ; 
PUT  EniT  fill  I  III  II  I  I  I  III  I  II  III  Mllll')  (COL  (91),  A(  30)  J? 

PUT    SKIP     (01 ; 

PUT    EDIT     ('XXXXXXXXX    XXXXXXXXX")     (COL    (101),     A(20)l: 
PUT    SKIP     (0); 

PUT    EDIT     (•  tDiUtlUHHVH  ■)     (COL(lll).     A(IO)); 
END; 
PUT    SKIP; 

PUT    EDIT     ((•='     00     I    =     10    TO    1201)      (COL(I),    A(ll); 
PUT    EDIT     ('FREQUENCY')     (C0L(2).     Al9))i 
PUT     FDIT    (FREQIK, 11)     (F(6) ) ; 
DC    N    =    ?.    TO    11; 

PUT    EDIT    (FREQIK.NM     (X(4),     F16I); 
END; 
PUT    PACE; 

end; 
end;      /*  regin  block  */ 

END  HAP2; 
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HAP3:  PROCEDURE  CPTICNS  (►'AIM); 


/#  ERP  _  M/iP3  PY  SHOMAS  M.  KUNTZ 

/*  OEPT.  OF  REGIONAL  £  COMMUNITY  PLANNING 

/«  KANSAS  STATE:  UNIVERSITY 

/a*** ***»******#******************************************** ********** 

/* 

MAPH  IS  A  SPECIALIZED   PROGRAM  OF  THE  ERP  PACKAGE  AND  IS 
DESIGNED  TO  MAP  VARIABLES  THAT  ARE  EXTRACTED  BY  OCCURENCE.   IT 
PRODUCES  A  SPATIALLY  RFLATEI).  TUG  DIMENSIONAL  MAP  INDICATING  THE 
CELLS  IN  WHICH  DATA  IS  PRESENT.   THE  USER  MAY  MAP  UP  TO  TEN 
ATTRIBUTES  OF  A  SINGLE  VARIABLE.   THE  FOLLOWING  CARDS  ARE 
REQUIREO  BY  THF  USER  TC  ESTABLISH  THE  BOUNDARIES  CF  THE  SITF. 
INAX  =  (NUMBER  OF  ROWS) ; 
JMAX  =  (NUMBER  OF  COLUMNS); 
HE  USER  MUST  ALSO  PROVIDE  THE  REFERENCE  NUMBER  OF  UPPER  LEFT 
CELL  (NOTC  -  THIS  CEIL  MLST  HAVE  THF  LOWEST  ROW  AND  COLUMN  VALUE). 
/*  THE  USER  MUST  ALSO  PROVIDE  THE  INCREMENTAL  DIFFERENCE   BETWEEN 

/*  GRID  LINES. 

I  AXIS  =  (ROW  NUMBER  OF  UPPER  LEFT  CELL); 
JAXIS  *  (COLUMN  NUMBER  OF  UPPER  LEFT  CELL); 
INCREMENT  =  (INCREMENTAL  DIFF.  BETWEEN  GRID  LINES); 


/* 

/*' 
/* 

/v 

/* 
/* 
/* 
/* 
/* 
/* 


/* 

/* 

/■■ 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 


ADDITIONAL  CARDS  FOR  PRINTING  OUT  APPROPRIATE  TITLES  AND  A 
LEGEND  MUST  ALSO  BE  SUPPLIED  BY  THE  USER.   THF  VAPIABLE  NAME 
•TITLE'  IS  A  CHARACTER  STRING  OF  50  CHARACTERS  AND  APPEARS  AT  THE 
TOP  OF  THE  MAP.   ITS  INTENDED  PURPOSE  IS  TO  DESIGNATE  THE  NAME  OF 
THE  VARIABLE  BFING  MAPPED  OR  THE  PROJECT  TITLE.   THE  VARIABLE 
NAME  'LABEL  IS  ALSO  A  CHARACTER  STRING  CONTAINING  A  MAXIMUM  OF  50 
CHARACTERS.   IT  APPFARS  DIRECTLY  BELOW  THE  MAP   AND  MAY  BE  USED 
TO  IDENTIFY  THE   DATA   BEING  MAPPED.   THE  NAME  'LECt*)'  IS  AN 
ARRAY  OF  TEN  CHARACTER  STRINGS  EACH  50  CHARACTERS  LONG.   EACH 
VARIABLE.   EACH  VARIABLE  BEING  MAPPED  IS  ASSIGNED  A  VALUE  BETWEEN 
1  AND  10  WHICH  CORRESPONDS  TO   THE  LEG  ELEMENT  (LEG(l),  LEGC2JJ. 
IF  LESS  THEN  10  ATTRIBUTES  ARE  USED  THEN  THE  UNUSED  ELEMENTS 
MUST  BE  A  BLANK  ENCLOSED  AS  SINGLE  QUOTATION  MARKS  SUCH  AS: 

LEG (10)  =  *  * : 
THE  VARIABLE  NAMF  <LINE(*)<  IS  A  5  ELFMENT  ARRAY  THAT  PERMITS 
THE  USER  TO  PROVIDF   UP  TO  5  LINES  OF  SPECIAL  INFORMATION  OR 

THE  COMPUTER  PRINTOUT.   IF  THE  USFR 
OF  THE  LINE  OPTION  HE  MUST  INSFRT  A 


/*  PROJECT  IDENTIFICATION  ON 
/*  DUES  NOT  WISH  TO  MAKE  USE 


/*  CARD  WITH  THE  FOLLOWING  STATEMENT: 


/* 
/* 

DCL 
DCL 
DC  I 


LINEd  )  =' 


(  IMAX,  JMAX)  FIXED  BIN; 

(AXIS,  JAXIS,  INCREMENT)  FIXED  BIN; 

(TITLE,  LABEL,  LEG(lOl)  CHAR150); 
DCL  LINEI5)  CHAR(IOO);  __ 

/***«.****#x.  ***********************************************  *********** 


/*  USER  CARCS  FOLLOW  */ 

IMAX  = 
JMAX  = 


IAXIS  * 
JAXIS  ' 


/*  USER  SUPPLIED  INFORMATION  */ 
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INCRf  FNT  = 
CELL_i\UMB  = 
TITLfc  = 
LABELS*  )  = 
LINE(*I  = 


BEGIN; 

OPEN  FILE  (SYSPRINT)  PAGES! Z E ( 08 ) ; 

Df.l  (X,  Y)  FIXED  BIN; 
IKL  ATTR  FIXED  BIN; 
OCL  REMAP!  IWAX.JMAX I  FIXEC  BIN  INITIAL  ( I IMAX* JMAX ) 0 ) ; 

DCL  (1ROW,  JCOL)  FIXED  BIN; 

OCL  (ROWIIMAX),  COLIJMAXI)  FIXED  BIN; 
OCL  FREQUOI  FIXED  BIN  I  N  IT  I  AL  (  (  10  )  0  )  ; 
ON  ENDFILE  (SYS1M)  GO  TC  NAP; 


SEARCI- 
SEEK 


HUNT 


FOUND 


IRCW  =  IAXIS; 

l  =  l; 

GET  EDIT  (X )  (COL  (1),  F(4) )  ; 

IF  X  =  IROW  THEN  GO  TO  HUNT; 
IROW  =  I ROH  +  INCREMENT; 

I  =  I  +  i; 

GO  TO  SFEK; 
JCOL  =  JAXIS; 

GET  FDIT(Y)  (CCLI6),  F(4) ) ; 
DO  J  =  1  TC  JMAX; 

IF  Y  -=  JCOL  THEN 

JCOL  =  JCOL  <■  INCREMENT 
ELSE  GO  TO  FOUND; 
END: 
GET  EDIT  (ATTR)  (F(2) ) ; 
REMAP1I,J)  =  ATTR: 
FREQ(ATTR)  =  FREO(ATTR) 
GO  TO  search; 
IROW  =  IAXIS; 
JCCL  =  JAXIS; 
I  =  1  TO  I  MAX; 

pcwui  =  irow: 

irow  =  ircw  ^   increment; 

END; 

J  =  1  TO  jmax; 

col(J)  =  jccl: 

JCOL  =  JCOL  ♦  increment; 

END; 


I 


MAP: 


00 


no 


/*  INFORMATION  IS  TKANSLATFD  INTO  CHARACTER  SETS  AND  OUTPUTTEO  IN 
/*  THE  MAP  FCRMAT. 
OCL  SYM(IO)  CHAR1 I  I  ; 

SYMU)  -  '.' 

SY"(?) =  •  I1 

SYM( 31= '-• 

SYM(<V)  =  M' 

SYMI 5)='X< 

SYM(6)=»0" 

SYM( 7)='=* : 
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SYM(R)  =•  I  • : 
sym( 9)  =  'x' ; 

SYM 10)='*'  ; 

/*  ARRAY!*, *l  IS  A  TEMPORARY  STORAGE  SPACE  WHICH  IS  USED  TO 
/*  THE  COLUMN  REFERENCE  NUMBERS  VERTICALLY. 

DCL  ARRAY14.JMAXI  FIXED  BIN  INITIAL  t  J  4*  Jf"AX)  0)  ! 
OCL  (ITEMP,  NFIX141)  FIXED  BIN! 
DO  N  =  1  TO  JMAX  BY  2; 
ITEMP  =  COLIN)  ! 
DC  I  *  I  TO  3! 

NFIX(I)  =  !TEMP/10**14-1  ); 

ITEMP  =  ITEMP-(NFIX( I  1*10**14-1 )) ; 

ARRAYd  ,N)  =  NFIX!  I  ); 

END; 

ARRAY(4,N)  =  ITEMP! 
END! 


PRINT 


/*  BEGIN  PRINTING  MAP 

PUT  EOI TtTITLE )  (A150) )  ! 
PUT  SKIP! 
OC  V    =  I  TO  4! 

PUT  SKIP  EDIT  ('  ■)  IC0H6),  A(l))! 
DO  N  =  1  TO  J"AX  BY  2! 
PUT  EDIT  !ARRAY(M,N))  (F(2)l! 
END! 
END! 

PUT  EDIT  ((■*'  DO  I  *  6  TO  JMAX  +  9)1  ICOL(I),  All))! 
PUT  EDIT  !'*',  ■*•)  (C0LI6),  All),  C0UJMAX*9J,  Ail)); 
DO  I  =  1  TC  IMAX! 

PUT  EDIT  (POW  (I),  •*')  IC0U1I.  FCI,  XII),  All))! 
DC  J  =  1  TC  JMAX! 

IF  REMAP! I, J)  >  5  THEN  DO! 

PUT  EDIT  (SYMI6))  (CCLU  +  7),  Ad)); 
END! 
DO  N  =  I  TO  5! 

IF  REMAP! I, J)  =  N  THEN  DO! 

PUT  EDIT  (SYMIN))  (COLIJt-7),  Add! 
END! 
END! 
END; 
PUT  SKIP(O) ; 
DC  K  =  7  TC  10; 

CO  J  =  1  TC  JNAX; 

IF    REMAPd  ,J)     >=    N    THEN    DC! 

PUT    EDIT     (SYM(NI)     (COl.U  +  7).     All)); 

END; 

END; 
PUT    SKIPIO) ; 

FND; 

PUT    EDIT     !'    *     ',     ROWdl)     (COL  (J  MAX+O)  ,    A{:)),    F14II; 

END;  I*     I     LOCP       */ 

PUT    FD1T    (•*<,     ■*•!     1CULI6I,    Ad),    C.f)L  (JMAX+9)  ,    Ad)): 
PUT    EDIT     ((•*'     DO     I    =    6    TO    JMAX    +    91)     (COLIII,     All)); 
DO    M    =     1    TO    4; 
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PUT     SKIP    EDIT     (<     "!     !COLtM,     A  t  1  J  )  ; 
DO    N    =    1    TO    JMAX    RV    2; 
PUT    rD IT    (ARRAY(K.N) )     (F(2!  !  5 
END; 
END; 
PUT     SKI?<2)     EDIT     (LA13LU     (A(50)l; 
PUT    SKIP(2) : 
IF    LINE(l)    -i=     •     '     THEN    DO; 

PUT    EDIT    (  I.INEUI  I     (AtlOO)  )  ; 
DO   N  *    2   TO   5! 

put  edit  (line(n)i  (coli50),  a  (  50  )  )  ; 

end; 
eno; 
put  skip(3>  edit  ('legend'!  !c0l(53)t  a150)); 

DO  N  =  I  TO  10; 

IF  LEGIN)  ■>=  '  *  THEN!  DO; 

PUT  EDIT  [N,  •  =  ',  LEGINJI  ICOLI50),  E(2I,  A(3),  A(50)l; 
END; 
END! 
PUT  SK[P(3)  EDIT  ('FRFCUENCY  DISTRIBUTION  OF  DATA  CELL  VALUES  IN  EAC 
H  LEVEL • )  (AllOOJ ) ; 

PUT  SKIPC2I  FOIT  ('LEVELS        0')  (X(l).  A(1<W); 
DO  NLEVEL  «  I  TO  1C; 

PUT  EDIT  (NLEVEL)  (X(8),  F(2)); 
END; 
PUT  SKIP(l); 
PUT  FOIT  <(•='  DO  I  =  10  TO  120))  (COL(I).  A(ll); 


00  N  M  TO  5; 
IF  N  =  3  THEN 

DO; 

PUT  SKIP  EDIT  ('SYMBOLS')  (Xll),  A(7)); 
END; 
ELSE  CO! 

PUT  SKIP; 
END; 

PUT  EDIT  (• I  II  III  II  I  +++++♦+++  XXXXXXXXX  OOOIK 

0000  COOCOOOCO  COCOCCOOO  COCOOOOOO  000000000')  (CCL(21),  A(IOO)); 
PUT  SKIP  (0); 
PUT  EDIT  (>=========  =========  =========  =========1)  (C0H81), 

At  Ml  )  ; 
PUT  SKIP  (0); 

PUT  EDIT  ('  II  I  Ml  III  III  I1MII  III  I  III  II')  tCOLt  91),  At  30)); 
PUT  SKIP  to); 

PUT  EDIT  t'XXXXXXXXX  XXXXXXXXX')  (COL  (101),  A(20)li 
PUT  SKIP  (0); 

PUT  EDIT  («##*«*####•)  (COLtllll,  A(10)l; 
ENCJ 
PUT  SKIP; 

PUT  EDIT  ((•---'  DO  I  =  10  TO  120))  (COLtl),  Atll); 
TOT  =  o; 
DO  N  =  1  TO  10; 

TOT  =  FRFQ(N)  +  ICT; 
END; 
NZ1LCH  =  <JMAX*!MAX)-TOT: 
PUT  EDIT  ('FREQUENCY',  NMLCH)  (C0I.I2),  A  (  9 ) ,  Ft6)); 
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nc  N  -  1  TO  10; 

CUT  EDIT  (FRECXNI  1  (X(4),  F(6)  I; 

END; 
END;     /*  END  BEGIN  BLOCK   */ 
END  MAP  3! 
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MAPV.    PROCEDURE    OPTIONS    ("MM; 

/*    ERP    -    MAP'.  BY    THOMAS    M .    KUNTZ  * ' 

.%  DEPT.    OF    REGIONAL    AND    COMMUNITY    PLANNING    */ 

,*  KANSAS    S1ATE    UNIVERSITY  */ 

/»***«**»*♦**+******#********************«****************************' 

•  */ 

/*     MAP'*  IS  A  SPECIALIZED  PROGRAM  OF  THE  ERP  PACKAGE.   IT  IS  */ 

/*  OFSIGNFD  TO  MAP  RESOURCE  DATA  EXTPACTED  ACCORDING  TO  THE  NUMBER  */ 

/*  OF  OCCUPFNCES  PER  CELL.   AS  MANY  AS  35  ATTRIBUTES  OF  A  VAPIABLF  *) 

I*    MAY  BE  CAPPED.   THE  PROGRAM  WILL  ONLY  ACCEPT  INTEGER  VALUES  */ 

/*  BETWEEN  I  AND  10.   THE  PROGRAM  WILL  PRODUCE  A  SET  CF  SPATIALLY  */ 

/*  RELATED,  TWO-DIMENSIONAL  MAPS  (ONE  FOR  EACH  ATTRIBUTE)  ILLUSTRAT-  */ 

/*  ING  THE  DENSITY  PATTFRN  OF  EACH  RESOURCE  AND  THE  FREQUENCY  *l 

/*  DISTRIBUTION  FOR  EACH  VALUE.   THE  FOLLOWING  CARDS  ARE  REQUIRED  */ 

/*  BY  THE  USER  TO  ESTABLISH  THE  BOUNDARIES  OF  THE  SITE.  */ 

/*            [MAX  =  (NUMBER  OF  ROWS);  *' 

/#            JMAX  =  (NUMBER  CF  COLUMNS);  *' 

/*            KMAX  =  (NUMBER  OF  ATTRIBUTES);  *■ 

/*  THE  USER  MUST  ALSO  PROVIDE  THE  RFFERENCE  NUMBER  OF  UPPER  LEFT  *, 
/*  CELL. (NOTE  -  THIS  CELL  MUST  HAVE  THE  LOWEST  ROW  AND  COLUMN  VALUE).*, 
/*  THE  USER  MUST  ALSO  PROVIDE  THE  INCREMENTAL  DIFFERENCE   BETWEEN 

/*  GRID  LINES.  ' 

/*  IAXIS  =  (ROW  NUMBER  OF  UPPER  LEFT  CELL);  *' 
/*            JAXIS  =  (COLUMN  NUMBER  OF  UPPER  LEFT  CELL); 

/*            INCREMENT  =  (INCREMENTAL  DIFF.  BETWEEN  GRID  LIMES);  *. 

/* 

/*       ADDITIONAL  CARDS  FOR  PRINTING  OUT  APPROPRIATE  LABELS  MUST 

/*  ALSO  BE  SUPPLIED  BY  THE  USER.   THE  VARIABLE  NAME  'TITLE'  IS  A      *. 

/*  CHARACTER  STRING  OF  SO  CHARACTERS  USED  TO  IDENTIFY  THE  VARIABLE    * 

/*  NAME.'  'LARFU*!'  IS  A  VARIABLE  NAME   DESCRIBING  THE  SPECIFIC 

/*  ATTRIBUTE  OF  THE  VARIABLE  BEING  MAPPED.   THERE  IS  ONE  LABEL  FOR    * 

/*  EACH  ATTRIBUTE.   THF  MAXIMUM  LENGTH  IS  100  CHARACTFRS. 

/*  THE  USER  CARD  SHOULD  FCLLOW  THE  FORMAT  BELOW:  * 

/*  TITLE  =  'NAME  OF  VARIABLE';  * 

/#  LABELt*)  =  'NAME  OF  ATTRIBUTE';  .* 

/*  THF  ASTERISK  (*)  REPRESENTS  THE  SEQUENTIAL  NUMBER  OF  EACH 

/*  *TTRTHE  LSFR  MAY  ALSO  USE  UP  TO  FIVE  LINES  ON  THF  COMPUTER         * 

/*  PRINTOUT  TO  PROVIDE  INFORMATION  PFRTINENT  TO  HIS  SPECIAL  PROJECT 

/*  OR  STUDY.   THE  FIRST  LINE  MAY  CONTAIN  11°  TO  100  CHARACTERS  AND 

/*  THE  REMAINING  FOUR  LINES  CONTAIN  50  CHARACTERS  EACH.   IF  THF  USER  * 

/*  DOES  NOT  WISH  TO  USE  TEE  LINE  OPTION  HE  MUST  USE  THE  STATEMENT:    * 

/*  LINE (I )  =  '     •;  * 

* 
/* 

/*****#******************»********************************** ********** 
DCL  UMAX,  JMAX,  KMAX)  FIXED  BIN; 
DCL  (IAXIS,  JAXIS,  INCREMENT)  FIXED  BIN; 
DCL  TITLE  CHARI50)! 
DCL  LABELM5I  CHARUOO): 
/♦,M  ************ ********■<'*************************** ********* 

/*  USER  CARDS  FOLLOW  */ 

I  MAX  = 
JMAX  = 
KMAX  = 


[AXIS  = 
JAXIS  = 
INCREMENT  = 
TITLE  ■ 
LABEL!*)  = 
LINE!*)  = 


/*  USER  SUPPLIED  INFORMATION  */ 


( ( IMAX*JMAX*KMA 


x )  o  I  : 


SEARCH: 
SEEK: 


HUNT: 


FOUND; 


-■    IAXIS; 
■  I; 

GET  EDIT  (X  )  (COL  (1  )  ,  F  (  <t )  )  ; 
IF  X  =  I  ROW  THEN  CO  TO  HUNT; 
IROW  =  IROW  *  INCREMENT; 

1=1+1; 

GO  TO  SEEK; 
JCOL  =  JAXIS; 

GET  EDITIYI  (CCL(6).  FK1  I  ; 
00  J  =  1  TO  JMAXi 

IF  Y  ->=  JCOL  THEN 

JCOL  =  JCOL  +  INCRFMENT ; 
ELSE  GO  TO  FOUND; 

END; 

DO  K  =  1  TO  KMAX; 

GET  EDIT  IATTRIK ) )  (F(2) I ; 


REMAP!  I,  J,K) 
FREOIATTR(K) ) 
END; 


=  ATTR(K) ; 

=  FREQ(ATTR  IK)  )+l: 


GC  TO  SEARCH; 
MAP:  IROW  ■  IAXI S; 
JCOL  =  JAXIS; 
DO  I  =  I  TC  rNAX; 

ROWt I  )  =  I  ROW; 

IROW  =  IP  OK1  + 

END; 
DO  J  =  1  TO  JMAX 

COL(J)  =  JCOL 

JCOL  =  JCOL  i- 

END; 


INCREMENT: 


IKCREfENTi 


/*  INFORMATION  IS  TRAKSLATED  INTO  CHARACTER  SETS  ANO  OUTPUTTED  IN 
/*  THE  MAP  FORMAT. 
OCL  SYMUOI  CFARd  )  ; 

SYMt  1)  =  '.' 

SYM2)  =  '  I  ' 

SYN(T)='-' 

SYMt  <,)=•  +  • 

SYM(5)='X« 

SYM(6) ='0" 
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SYM(  71  =  •=• ; 
SYM(B)=" I  * ; 
SYf  C3)  =  '  X'  ; 
SYM(  10)  =  ,«' ; 
DCL  ARRAYK,J!*AX>  FIXED  BIN  INITIAL  t  ( <<*  JMA  X  )  0  )  ' 
DCL  IITEMP,  NFIXK)  FIXED  BIN! 
DO  N  =  1  TO  J1AX  BY  2; 
ITEN'P  =  COLIN)  ; 


DC  I  =  1.  TO  3; 

NF1XI  I  )     =    ITFMP/10**(<t- 
ITEMP    ■     ITEMP-lNEIXl  I  )■ 
APRAYU  ,N!     =    NFIX!  I  ); 
END; 

ARRAY(4,NI     =     ITEMP; 

E  ND  ; 


1 1 ; 

10**1  4- 


i ) ) 


/*  BEGIN  PRINTING  MAP 
DO  K  =  1  TO  KMAX; 

PUT  ECITIT  ITLE)  [ A ( 50 » ) i 

PUT    SKIP; 

DO    M    =    1    TC    4; 

PUT    SKIP    EDIT     (•     •)     (COL  (6).     All)); 
DO    N    =     1    TO    JMAX    BY    2; 
PUT    EDIT     (  ARRAYC.NI  )     (F(2))! 

END; 
END; 

PUT  EDIT  ((•*'  DO  I  =  6  TO  JI"AX  +  91)  (CUL(I),  A  (  1  I  I  ; 
PUT  EDIT  («*«,  ■*')  (CCLlO),  All),  C0L{JMAX*9)i  A(ll); 
fin  I  =  1  TC  SMAXi 

PUT  EDIT  (ROW  (I),  '*')  (CCLll),  F(4»i  Kllli  A (II) 5 
DO  J  --  1  TO  JMAX; 

IF  REMAPl I ,J, K)  >  5  THEN  DO  5 

PUT  EDIT  (SYM(61)  (COLU  +  7),  All)): 

end; 

DO  N  =  1  TO  5; 

IF  REMAPU  ,J,K)  =  N  THEN  DO; 

PUT  EDIT  (SYMIN))  (COLU+7),  A(ll); 
END; 

end; 

END; 
PIT  SKIP(O); 
00  N  =  7  TG  10; 

DC  J  =  I  TO  jr'AX; 

IF  RFMAPl I ,J,K)  >=  N  THEN  DO; 

PUT  EDIT  (SYM(N)I  (CnUJ+7),  All)); 

end; 
END; 
PUT    SKIP(O) ; 
END; 
PUT    EDIT     ('     *     ',     ROW(I))     (COL(  JMAX  +  fl)  ,    A(3),    FI4I): 
END;  /*    I    LOCP      */ 

PUT    EDIT    («*',     '**)     (COLif.l,     All),     CUL  [J  VtyX  +  9 )  ,    AID); 
PUT    EDIT     (('*«     DO     I    =    6    TO    JMAX    <-    9  >  )     (COM  lit    A(ll); 
DC    M    =    I    TO    '<; 

PUT     SKIP    EDIT     (<     •)     IC0L16),     All)); 
DO    N    -    1     TC    JCAX    BY    2; 
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PIT  EDIT  ( ARRAY! f  ,M I  (F(2I I ; 
END; 
ENC; 
END; 
PUT  SKIP12)  EDIT  (LARFL(K))  (A(50))i 
PUT  SKIPt? ) ; 
IF  LINE(l)  -■=  •  '  THEN  DOi 

PUT  EC  IT  (L  INCl  1)  )  14(  100  11  ; 
DO  N  =  2  TC  5  ; 

PUT  FDIT  (LINEIN1)  ICCL150),  A(50)); 
END; 
END; 
PUT  SKIPI3)  EDIT  ('FREQUENCY  DISTRIBUTION  nF  DATA  CFUL  VALUES  IN  FA 
H  LEVEL'  )  ( At  100)  ); 

PUT  SKIPI2)  EDIT  ('LEVELS        0')  (Xtl),  MHI); 
DO  NLEVFL  =  1  TO  1C; 

PUT  EDIT  (NLEVEL)  (XI8),  F(2)  )  ; 
END; 
PUT  SKIP (I  I  ; 
IF  X  =  3  THEN 
CO; 

PUT  SKIP  EniT  ('SYMBOLS')  (Xtl),  A(7>); 
END; 
ELSE  co; 

PUT  SKIP; 

PUT  EDIT  (• 1111111  II  tuttt  +  tt  XXXXXXXXX  0000' 

0000  000000000  OOOOCCOOO  CCCCCCOOO  0000000)0')  (C0L121),  A(IOO)); 
PUT  SKIP  (0) ; 
PUT  EDIT  (•=========  =========  =========  =========<)  (C0L!81), 

A(40) ); 
PUT  SKIP  10) ; 

PUT  EDIT  Clllllllll  lllllllll  lllllllll')  (C0L(91),  AOOIIl 
PUT  SKIP  (0| ; 

PUT  EDIT  ('XXXXXXXXX  XXXXXXXXX")  (COL  (101),  A(20)); 
PUT  SKIP  (0) ; 
PUT  EDIT  (•##*«#«»#«»)  (COL(lll),  AI101); 

end; 

PUT  SKIP; 

PUT  EDIT  ((•='  DO  I  =  10  TC  120))  (COL(l),  A(l)); 

TOT  =  o; 

DO  N  =  1  TC  10; 

TOT  =  FRFOtl  )  +  TCT; 

END; 
NZILCH  =  ( JWAX*I^AX)-TCT; 

PUT  EDIT  (FREQINI)  ( X ( 4 ) ,  FI61); 

FND; 
PUT  PAGE; 
END; 
END  «AP4; 
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The  application  of  a  computerized  geo-information  system  to  environ- 
mental land  use  planning  has  several  advantages.  First,  it  consolidates 
information  from  many  sources  Into  a  single  spatially-related  data  base. 
This  means  the  planner,  or  any  potential  user  of  the  information,  does 
not  have  to  seek  out  data  from  many  unrelated  sources  as  has  been  the 
traditional  approach.  Secondly,  the  data  base  which  has  been  extracted 
from  maps  and  aerial  photographs  of  varying  scales  1s  displayed  in  a  uni- 
form format,  thus  enabling  the  user  to  analyze  the  data  more  readily. 
Finally,  it  enables  the  user  to  analyze  data  on  a  regional  scale  with  a 
level  of  discreteness  that  is  virtually  impossible  with  traditional  tech- 
niques involving  overlays. 

The  computer  printout  is  a  graphic  display  of  the  data  in  map  form. 
Each  cell  is  represented  by  a  printed  character  that  is  1/8"  x  1/10", 
thus  causing  an  elongation  in  the  vertical  dimension.  The  user  must 
remember  that  the  data  has  been  extracted  on  a  square-cell  basis  and  that 
the  distortion  is  only  graphic.  Further  distortion  may  occur  from  using 
the  cell  system  to  store  linear  elements  such  as  highways.  Diagonal 
elements  will  intersect  adjacent  cells  in  such  a  way  that  they  will 
appear  to  deviate  from  a  straight  line  on  the  printout.  Additional  dis- 
figuration will  occur  as  a  result  of  the  UTM  grid  deviating  from  true 
north.  The  problem  is  related  to  the  selection  of  cell  size.  As  cell 
size  decreases,  a  linear  feature  will  more  closely  approximate  its 
true  configuration. 

Adjacent  to  the  border  of  the  study  area  are  the  four-digit  cell 
identification  numbers.  A  user  interested  in  the  data  of  a  particular 
cell  or  cluster  of  cells  may  readily  locate  the  site  by  using  a  U.S.G.S. 
map  keped  with  the  UTM  grid  lines. 
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Appearing  at  the  top  of  the  printout  is  the  name  of  the  variable  or 

component  of  the  resource  Information.  Below  the  map  the  attribute  or 
specific  data  description  1s  identified  and  is  followed  by  the  study 
title.  Below  the  title  is  a  legend  which  denotes  what  each  particular 
symbol  represents.  There  are  two  general  categories.  In  the  first  case 
a  symbol  may  represent  an  attribute  and  simply  indicates  that  the  resource 
is  present  within  the  cell.  For  example,  the  map  titled  "Surface  Drainage 
Patterns,"  (page  62)  intermitent  streams  are  represented  by  symbol  1,  (a 
period)  whereas  a  major  river  1s  represented  by  symbol  9,  (an  over 
printing  of  an  0,  =,  1,  and  x).  In  the  second  case,  a  symbol  may  repre- 
sent a  value  range  for  an  attribute  as  1s  the  situation  where  data  is 
extracted  by  percent  of  cell.  The  value  ranges  that  the  symbols  repre- 
sent can  be  found  below  the  line  "Value  Range  of  Each  Level."  Topographic 
Relief  (page  50)  is  mapped  in  10  levels  between  the  lowest  centroid  ele- 
vation and  highest  centroid  elevation.  Data  extracted  by  percent  of  cell 
is  mapped  in  10  levels  between  0  and  100  percent.  Directly  above  the 
symbol  block  is  the  percentage  of  the  total  site  occupied  by  the  resource. 
The  last  line  on  the  printout,  "Frequency,"  is  a  record  of  the  number  of 
occurrences  of  each  symbol. 
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FURTHER  RESEARCH 


The  geo-1nformat1on  system  developed  in  this  study  is  merely  an 
accounting  procedure.  Information  is  extracted  from  many  sources; 
it  is  stored  in  a  computer  in  a  spatially  relatable  manner,  and  then 
may  be  retrieved  in  the  form  of  a  two-dimensional  graphic  display. 
Although  the  system  is  of  great  benefit  to  the  land  use  planner  due 
to  the  ease  of  which  Information  may  be  accessed  the  system  has  much 
greater  potential . 

The  analysis  of  data  contained  in  the  information  system  through  a 
decision  making  model  is  the  next  logical  phase  in  the  development  of  a 
computerized  geo-information  data  bank.  Such  a  model  would  involve  the 
development  of  a  criterion  for  optimum  land  use  such  as  proposed  by  Ian 
McHarg  (HcHarg,  Twin  Cities  Study,  196S). 


FIGURE  6  DECISION  MODEL 


The  computer  procedure  would  search  each  cell  for  the  appropriate 
criterion  and  then  provide  a  graphic  display  ranking  each  cell's  suita- 
bility for  the  intended  land  use.  The  development  of  such  a  model  would 
be  appropriate  for  further  application  of  this  study. 
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Another  area  for  further  study  is  the  matter  of  selecting  cell 
size.  Currently  the  method  of  selecting  an  appropriate  cell  size  may 
be  at  best  described  as  arbitrary.  It  would  seem  that  a  more  scientific 
system  could  be  developed.  Perhaps  a  method  of  sampling  various  cell 
sizes  and  performing  a  statistical  test  for  significant  difference 
could  be  developed  as  a  supplemental  research  project. 
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The  Computer  System 

The  modern  computer  is  cm  electronic  device  capable  of  mani- 
pulating enormous  quantities  of  data  and  completing  numerous  mathe- 
matical calculations  of  extreme  complexity  at  tremendous  speed.  The 
computer  system  consists  of  five  basic  components.  The  Central 
Processing  Unit  performs  the  basic  mathematical  operations  and  controls 
the  interaction  of  the  other  elements  in  the  system.  The  Compiler  is 
essentially  a  device  for  facilitating  communication  between  man  and 
machine.  It  consists  of  a  set  of  instructions  for  translating  a 
Program  Language,  such  as  Fortran,  PLI,  COBAL,  etc.,  into  Machine 
Language.  The  Memory  is  the  limiting  factor  in  determining  how  large 
a  problem  the  computer  can  handle  at  one  time.  The  Input  and  Output 
devices  (often  referred  to  as  1/0)  make  convenient  provision  for 
reading  data  and  instructions  into  the  computer  and  for  getting  the 
analysis  out.  A  number  of  different  types  of  input  and  output  equip- 
ment are  available.  For  simplicity,  the  geo-information  system  is 
designed  to  input  data  on  the  Card  Reader  and  output  information  on  the 
Line  Printer.  The  card  reader  is  the  most  common  input  device.  There 
are  however,  numerous  advantages  of  using  other  input  devices  such  as 
Tape  or  Disk.  The  user  may  obtain  assistance  from  the  consulting  services 
at  the  computer  center  for  using  either  of  these  two  storage  modes. 

The  Punch  Card 

The  punch  card  has  spaces  for  small  rectangular  holes  organized 
1n  80  columns  and  12  rows.  These  holes  are  punched  on  a  Key  Punch 
Machine  which  resembles  an  electric  typewriter  and  which  may  be  found 
in  any  computer  center.  Information  is  recorded  from  right  to  left 
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FIGURE  7     BLOCK  DIAGRAM  OF  THE  COMPUTER  SYSTEM 
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FIGURE  8  THE  IBM  PUNCH  CARD 


90 
across  the  card  in  the  same  fashion  as  information  is  typewritten  across 
a  page.  The  80  columns  correspond  to  80  spaces  which  80  numbers, 
letters,  or  other  symbols  may  be  recorded.  Each  column  has  spaces  for 
as  many  as  12  punches,  one  in  each  of  the  12  rows.  The  rows  are 
numbered  0-9  with  the  9- row  at  the  bottom  of  the  card,  and  there  are 
two  additional  rows  (sometimes  called  the  11  and  12-row)  at  the  top  of 
the  card.  Numbers  are  recorded  by  punching  a  single  hole  1n  one  of  the 
0-9  rows  of  a  given  column.  Letters  and  punctuation  symbols  are  recorded 
by  punching  two  or  more  holes  in  a  given  column.  For  example,  the 
letter  A  consists  of  a  punch  in  the  very  top  row  (the  12-row)  and  1n  the 
1-row,  the  letter  B  of  a  punch  in  the  very  top  row  and  in  the  2-row,  etc. 
The  Card  Reader  translates  these  punches  into  binary  coded  information 
reading  each  column  at  a  time  and  checking  each  row  to  see  if 
there  is  a  punch. 

When  numerical  data  are  restricted  to  the  positive  integers  0-9, 
any  number  may  be  recorded  in  a  single  column  on  the  card,  and  the  card 
may  contain  as  many  as  80  different  pieces  of  numerical  information. 
However,  many  of  the  numbers  we  will  wish  to  record  on  punch  cards  will 
be  larger  than  one  or  may  have  decimal  values.  These  numbers  are 
recorded  by  grouping  columns  together  to  form  fields  of  the  desired 
width.  Thus,  a  punch  card  used  for  numerical  data  may  contain  80  fields 
of  one-column  width,  40  fields  of  two-column  width,  16  fields  of  five- 
column  width,  or  any  combination  of  fields  of  varying  width  as  long  as 
the  total  number  does  not  exceed  80.  The  field  width  to  be  used  in 
the  geo-information  data  bank  have  been  predetermined  and  are  speci- 
fied in  the  programming  section. 
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Organizing  The  Program  For  Computer  Analysis 

Each  program  consists  nf  a  set  of  program  decks  which  must  be 
in  the  proper  sequence  for  computer  analysis.  The  first  deck  is  a 
series  of  Job  Control  Language  (OCL)  cards  which  contain  information 
relating  to  the  execution  of  the  program.  The  second  deck  contains  a 
set  of  user  cards  in  which  the  user  defines  the  boundaries  of  the 
study  area  and  provides  certain  additional  information  required  for 
the  execution  of  the  program.  The  third  set  of  cards  is  the  program 
deck  which  contains  the  execution  cards  and  causes  the  program  to 
process  and  display  information  contained  in  the  data  bank.  This 
deck  is  proceeded  by  another  OCL  card.  The  fourth  series  of  cards  in 
the  data  deck  which  contains  the  extracted  information  relating  to  a 
particular  variable.  It  is  followed  by  a  final  OCL  card. 


DATA  DECK 


DCL 

ROGRAM  DECK 


"USER  DECK 


FIGURE  9  ORGANIZATION  OF  A  PROGRAM 
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Job  Control  Language 

The  Job  Control  Language  consists  of  a  set  of  control  cards  used 
to  direct  the  execution  of  the  computer  program  and  to  describe  the 
input/output  devices  used.     JCL  varies  with  different  systems  and 
computer  centers.     The  following  is  a  description  of  the  OCL  state- 
ments used  for  this  program.     The  user  may  wish  to  consult  a  JCL 
reference  guide  at  this  particular  computer  center  for  any  variations 
1n  format. 

1.  The  JOB  statement.     This  is  the  first  card  and  marks  the 
beginning  of  a  job. 

//SITEMAP        JOB  (        ssi\      .  acct  f  ,3,6,  ,1802,1 ,  ,0), 'name' 

The  above  JOB  statement  will  provide  for  an  output  of 
6,000  lines  printed  on  white  paper  at  8  lines  per  inch. 

2.  The  EXEC  (Execute)  statement.     This  card  follows  the 
JOB  card  and  names  the  program  or  procedure  to  execute. 

//ERP  EXEC  PL1LFCLG,PARM,PL1L='LINECNT=84' 
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3.  The  DD  (Data  Definition)  statement.  This  card  describes  each  data 
set  (in  this  case  a  deck  of  cards)  and  requests  the  allocation  of 
I/O  devices.  This  program  requires  two  DD  statements.  The  first 
card 

//PL1.SYSIN  DD  * 

follows  the  EXEC  statement  and  the  second  card 

//GO.SYSIN  DD  * 

precedes  the  data  deck. 

4.  The  Delimiter  (/*)  statement.  This  is  the  "end-of-file"  card 
for  marking  the  end  of  the  data  deck. 


/* 


It  should  be  noted  that  all  JCl.  cards  (except  the  /*  card) 
begin  with  a  //  in  columns  1  and  2,  followed  by  a  name  field. 
The  words  JOB  and  EXEC  begin  in  column  12. 
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The  User  Deck 
As  previously  stated  the  function  of  the  user  deck  1s  to 
define  the  boundaries  of  the  site  and  provide  additional  information 
necessary  before  the  execution  of  the  program  can  begin.  The 
Information  required  for  each  MAP  program  will  not  always  be  the 
same,  however,  the  variable  names  used  1n  each  program  are  constant. 
The  following  is  a  list  of  the  variable  names  and  a  description  of  the 
type  of  information  required. 


Variable  Name 
IMAX 
OMAX 
KMAX 

2MAX(*) 
CELLJUHB 
3IAXIS 
3JAXIS 

Increment 

TITLE 
LABEL(*) 

LINE(*) 

2LEG(*) 


Description 

Integer  value,  maximum  number  of  rows. 

Integer  value,  maximum  number  of  columns. 

Integer  value,  maximum  number  of  attributes. 

Integer,  array,  minimum  value  of  each 
attribute. 

Integer,  array,  maximum  value  of  each 
attribute. 

Integer,  total  number  of  cells  in  study 
area. 

Integer,  identification  number  of  upper 
most  row. 

Integer,  identification  number  of  lefthand 
col unm . 

Integer,  values  by  which  identification  number 
are  incremented. 

Character-string,  length  50,  name  of  variable. 

Character-string,  length  100,  name  of  each 
attribute. 

Character-string,  array,  length  100,  project 
description,  maximum  of  5  lines  of  print. 

Character-string,  array,  length  50,  name  of 
attribute  corresponding  to  coded  value,  maxi- 
mum 10. 
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1MAP1  only 
2MAP3  only 


3MAP3  and  MAP  4  only 

The  program  language  used  in  PL1.  It  1s  a  versatile  language 
and  1s  well  suited  to  the  type  of  operations  performed.  For  those  not 
familiar  with  the  programing  language  the  following  rules  are  provided 
to  aid  in  programing  the  user  deck. 

1.  All  statements  in  the  user  deck  should  begin  in  column  2. 
Do  not  use  column  1. 

2.  Do  not  go  past  column  72.  A  statement  may  be  continued 
on  a  second  card  beginning  in  column  2. 

3.  Each  statement  must  end  with  a  semi-colon  (;)• 

4.  Character-strings  must  be  enclosed  with  single  quotation 
marks  ('character-string'). 

5.  A  statement  beginning  with  a  /*  and  ending  with  an  */  is 
a  comment  card.  This  is  a  nonexecutable  statement  and 
is  used  to  aid  in  documenting  the  program. 

The  program  documentation  for  each  MAP  program  specifies  the 
required  user  cards.  The  user's  attention  is  called  to  the  statement. 

/*  USER  CARDS  FOLLOW  */ 

This  1s  the  last  statement  of  the  program  documentalon  and  is  proceeded 
by  the  user  cards. 
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The  following  examples  are  provided  to  Illustrate  how  the  user 
cards  for  each  MAP  program  arc  to  be  punched.  Assume  for  each  example 
that  an  Irregular  shaped  site  is  being  mapped.  The  maximum  number  of 
cells  in  a  column  (IMAX)  is  100  and  the  maximum  number  of  cells  in  a 
row  (JMAX)  is  75  and  the  total  number  of  cells  (CELLJUMB)  is  5350. 
The  upper  right  hand  cell  identification  number  is  100J00  and  the 
values  are  incremented  by  10. 

Example  1.  The  Use  of  MAPI. 

Let  us  assume  we  are  going  to  use  the  MAPI  program  to  record  the 
resource  variable;  topographic  relief.  Three  attributes  have  been 
extracted  from  each  cell:  centroid  elevation,  lowest  elevation,  and 
highest  elevation.  The  table  below  indicates  the  minimum  and  maximum 
values  of  each  attribute  over  the  entire  site. 

Attribute  Minimum  Value  Maximum  Value 

(1)  Centroid  Elevation         900'  1200' 

(2)  Lowest  Elevation  860'  1170' 

(3)  Highest  Elevation  930'  1220' 

The  program  requires  the  following  user  cards  to  be  provided. 

IMAX  =  100; 

JMAX  =  75; 

KMAX  =  3; 

CELL  NUMB  =  5350 

MIN(T)  =  900; 

MIN  2  =  860; 

MIN(3)  =  930; 

MAX(i  =  1200; 

MAX(2)  =  1170; 

MAX(3)  =  1220; 

TITLE  =  'TOPOGRAPHIC  RELIEF'; 

LABEL(l)  =  'ELEVATION  MAP  RASED  ON  CENTROID  ELEV. ' ; 

LABEL(2)  =  'ELEVATION  MAP  BASED  ON  LOWEST  ELEV.'; 

LABEL(3)  *  'ELEVATION  MAP  BASED  ON  HIGHEST  ELEV.'; 
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Additionally,  the  user  is  allocated  five  lines  of  print  to 
provide  information  relating  to  the  project.  Information  such  as 
project  title,  student  name,  course  name  typify  information  that 
may  be  used.  The  following  is  an  example  of  the  use  of  this  option. 

LINE(l)  =  'ENVIRONMENTAL  LAND  USE  PLANNING  STUDY'; 

LINE(2)  =  'JOHN  E.  DOE'; 

L1NE(3)  =  'DEPT.  OF  LANDSCAPE  ARCHITECTURE'; 

LINE  4  =  'KANSAS  STATE  UNIVERSITY'; 

LINE(S)  =  '  '; 

Note  that  if  all  5  lines  are  not  used  a  statement  must  be  included 
with  single  quote  marks  around  a  blank  column.  If  the  user  does  not 
wish  to  use  the  LINE  option  he  need  only  insert  the  statement; 

LINE(*)  =  '  '; 

which  will  cause  the  option  to  default. 

The  data  deck  contains  the  extracted  information  pertaining  to  a 
particular  resource  variable.  Each  card  is  a  record  containing  a  cell 
identification  number  followed  by  one  or  more  fields  of  data.  These 
fields  contain  the  resource  information  extracted  from  the  cell 
identified. 

The  MAPI  program  uses  a  5-column  field  width,  beginning  with 
column  11,  for  recording  each  attribute.  The  fields  are  numbered 
successively  and  corresponds  with  the  label  number.  As  many  as  11 
attributes  may  be  recorded  on  a  single  data  card.  A  typical  data 
card  for  data  extracted  for  Example  1  would  be  punched  as: 


m 


col  umn 

number 


1111111111 222222222233333333334 
1 23456789ai2345678S01 2345678901 234567890 


1025  4200 


1130 


1110: 


117CJ 


This  record  tells  us  that  1n  cell  number  1025  4200  the  centroid 
elevation  is  1130  feet,  the  lowest  elevation  is  1110  feet,  and  the 
highest  elevation  is  1170  feet.  MAPI  requires  a  data  card  for  each 
cell.  If  the  cell  is  not  included  in  the  study  area  the  fields  for 
the  attributes  are  left  blank. 


Example  2.  The  Use  of  MAP 2. 

The  MAP2  program  is  utilized  with  data  extracted  as  a  percent 

of  cell.  In  this  example  let  use  assume  we  are  mapping  soil  types  as 

a  resource  variable.  If  five  soil  classes  have  been  identified  and 

extracted  the  user  cards  required  are: 

IMAX  =  100; 

JMAX  =  75; 

KMAX  =  5; 

CELL  NUMB  =  5350 

TITLE  =  'SOIL  STRUCTURE  MAP'; 

LABEL(l)  =  'CLASS  1  SOILS'; 

LABEL  2  =  'CLASS  2  SOILS': 

LABEL(3)  =  'CLASS  3  SOILS': 

LABEL(4)  =  'CLASS  4  SOILS' 

LABEL(5)  =  'CLASS  5  SOILS', 

LINE(l)  -  'ENVIRONMENTAL  RESOURCE  PLANNING  STUDY'; 

LINE(2)  =  'JOHN  E.  DOE'; 

LINE(3)  =  'DEPARTMENT  OF  LANDSCAPE  ARCHITECTURE'; 

LINE(4  =  'KANSAS  STATE  UNIVERSITY'; 

LINE(5)  =  '   '; 
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The  date,  cards  follow  the  same  format  as  the  MAPI  program.  The 
first  9  columns  contain  the  cell  identification  number  followed  by  as 
many  as  14  fields  of  5-column  width  beginning  in  column  11.  Each  field 
contains  the  percent  of  coverage  for  each  attribute  in  that  cell.  The 
field  number  corresponds  to  the  LABEL  number  of  each  attribute.  If  the 
value  is  punched  in  the  last  two  columns  of  each  field  it  is  not 
necessary  to  punch  the  decimal  point.  If  no  value  is  recorded  in  a 
specified  field  then  the  value  is  assumed  to  be  zero.  A  typical  data 
card  may  be  punched  in  the  following  manner. 


col umn 
number 


1 
1234567890 


1025  4200 


11111111222222222233333333334 
2345578901 2345678901 234567890 


39 


45; 


The  above  card  indicates  that  252!  of  cell  number  1025  4200  is 
covered  by  class  1  soil,  0%   by  class  2  soil.  35%  by  class  3  soil, 
45%  by  class  4  soil,  and  0%  by  class  5  soil. 

Example  3.  The  Use  of  MAP3. 

The  HAP 3  program  produces  a  composite  map  of  attributes  extracted 
on  the  basis  of  predominate  type  or  may  be  used  to  map  data  on  the  basis 
of  occurrence  within  a  cell.  To  illustrate  the  use  of  the  program 
assume  that  the  soil  type  data  had  been  extracted  by  predominate  type 
rather  than  percent  of  cell  as  in  Example  2.  The  user  cards  required 
for  this  program  will  be: 


100 


IMAX  =  100; 

JMAX  *   75; 

IAXIS  ■  100; 

OAXIS  =  100; 

INCR  =10; 

TITLE  '   ' PEDLOGIC  RESOUP.CER'; 

LABEL  =  'SOIL  STRUCTURE  MAP';. 

LEG(l)  =  '  '; 

LEG(2)  =  'CLASS  1  SOILS': 

LEG(3)  =  '   '; 

LEG  4  =  'CLASS  2  SOILS' 

LEG(5)  ■=  '  '; 

LEG(6  =  'CLASS  3  SOILS': 

LEG(7)  -  '  '; 

LEG(8  =  'CLASS  4  SOILS'; 

LEG(9)  =  '  '; 

LEG(IO)  =  'CLASS  5  SOILS'; 

LINE(l)  =  'ENVIRONMENTAL  LAND  USE  PLANNING  STUDY'; 

LINE  2  =  'JOHN  E.  DOE'; 

LINE(3  =  'DEPARTMENT  OF  LANDSCAPE  ARCH.'; 

LINE(4)  ■  'KANSAS  STATE  UNIVERSITY'; 

LINE(5)  =  '   '; 

Note  that  there  is  only  one  LABEL  card  for  this  program.  The 
variable  name  'LEG'  is  used  to  print  out  the  appropriate  legend.  Each 
attribute  that  is  being  extracted  is  assigned  a  two-digit  integer  value 
that  corresponds  to  its  position  in  the  LEG  array.  In  this  instance 
class  1  soils  are  coded  as  02,  class  2  soils  as  04,  class  3  soils  as 
06,  class  4  soils  are  08,  and  class  5  soils  as  10.  The  program  is 
limited  to  mapping  ten  attributes  for  each  resource  variable.  Each 
unused  element  of  the  LEG  array  is  represented  by  a  blank  column 
between  single  quote  marks  ('  ). 

The  data  cards  contain  only  the  cell  identification  number  and 
the  two  digit  value  of  the  dominate  attribute  or  resource  found  in 
the  cell.  This  value  is  punched  1n  columns  10  and  11.  A  data  card 
for  the  cell  used  in  Example  2  would  be  punched  as  follows: 


101 


column 
number 


1111111111222222222233333333334 
1 234  5678901,2345678901 2345678901 234567890 


1025  4200 


The  predominant  soil  type  is  the  class  4  soil  which  corresponds  to 
the  eighth  position  of  the  LEG  array  and  is  therefore  coded  8.  Unlike 
MAPI  and  RAP2,  MAP3  requires  a  record  of  only  those  cells  from  which 
data  1s  extracted  or  in  which  the  resource  being  mapped  occurs.  This 
feature  facilitiates  extracting  data  by  occurrence  of  resources  which 
may  occupy  relatively  few  cells. 

Example  4.  The  Use  of  MAP4. 

The  MAP4  program  is  used  for  data  extracted  according  to  number 
per  cell.  It  functions  similar  to  the  MAPI  program,  however,  it  does 
not  require  a  data  card  for  each  cell.  The  program  assumes  that  if  a 
record  is  not  provided  for  a  cell  the  resource  does  not  occur  at  that 
cell  location.  The  program  is  also  limited  to  recording  integer  values 
from  1  to  10.  To  illustrate  the  use  of  the  program  assume  that  data 
relating  to  transportation  systems  are  to  be  mapped  according  to 
number  per  cell  and  the  following  attributes  have  been  extracted: 

1.  4- Lane  Divided  Highway 

2.  2- Lane  State  Highway 

3.  County  Road 

4.  Local  Serving  Street 

5.  Railroad 
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The  following  user  cards  are  required. 

IHAX  =  100; 
OMAX  «  75; 
KMAX  =  5; 

IAXIS  -  100; 

OAXIS  =  100; 

INCREMENT  -  10; 

TITLE  -  'TRANSPORTATION  ANALYSIS'; 

LABEL(l)  =  '4-LANE  DIVIDED  HIGHWAY'; 

LABEL(Z)  =  '2-LANE  STATE  HIGHWAY'; 

LABEL(3)  =  'COUNTY  ROAD'; 

LABEL (4)  =  'LOCAL  SERVING  STREET'; 

LABEL(5)  =  'RAILROAD'; 

The  program  will  produce  a  set  of  computer  maps,  one  map  for  each 
attribute,  which  graphically  depicts  the  density  of  each  transportation 
system. 

The  data  record  consists  of  the  cell  identification  number 
followed  by  as  many  as  35  2-column  fields  for  recording  the  number  of 
occurrences  of  each  attribute  within  the  cell.  The  field  number  of 
each  attribute  corresponds  to  its  respective  LABEL  number.  From  the 
above  example  a  cell  containing  a  4-Lane  divided  highway  and  2  county 
roads  would  be  punched  as  follows: 


column  1111111111222222222233333333334 

number     1 23456789pi  23,45,67,8°pi 2345678901 234567890 
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The  following  presenty  a  summary  of  one  classifica- 
tion, and  a  listing  of  data  as  mapped  on  individual  map 
sets,  1  liu  mapping  codes  for  each  item  in  itie  classifica- 
tion are  shown.  Area  data  arc  mapped  and  recorded  as 
areas;  only  the  point  locations  of  other  data  are  mapped. 
They  arc  recorded  as  mileages,  type1.,  or  numbers  of 
items.  , 

Land  Use  Classification 

(SUMMARY  SHEET) 
Agriculture 

Area:;: 

At— Orchard,  vineyard,  horticulture 

Ac— Cropland  and  pasture 
Point  data: 

Active  farms  (f) :  number 

Greenhouses  (g):  number 
Forest  Land 
Areas : 

Fb-Brushland 

Fn— Forest  stands 

Fp-Plantations  (10  acres  +) 
Water  Resources 
Areas: 

Wh— Hudson  River 

Ws— Other  rivers  and  streams  (100'  4 ) 

Wn— Natural  ponds  and  lakes  (1  acre  r) 

Wc~ Artificial  ponds  and  reservoirs  (1  acre  -f) 
Point  data: 

Natural  ponds  and  lakes  (n):  number 

Artificial  ponds  and  reservoirs  (c):  number 

Ponds  less  than  1  acre  in  size  (p):  number 

Take  and  reservoir  shoreline  (1):  miles 

Streams  (s):  miles 

Hudson  River  shoreline  (m):  miles 
Wetlands 
Areas: 

Wb— Bogs,  marshes 

Wbs— Salt  marshes 

Ww— Wooded  swamp 
Islands 
Areas : 

5  -Island  designation  (Used  as  prefix  for  land 
uses  occurring,  on  an  island 
Residential  Land  Use 
Areas: 

Rh-  -High  density  (50'  frontage  maximum) 

Rm -Medium  density  (50'-]  00'  frontage) 

Rl-Lo\v  density  (100'  !-  frontage) 

Re— Residential  estates  (larger  than  5  acres) 

Rr-  Rural  hamlet 
Point  data.: 

1  ligtwise apartment  buildings  (2):  number 

Trailer  parks  (v) :  number 
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Commercial  Land  Area 
Areas: 

Cu-Ccntrnl  business  district 

Cc— Shopping  center 

Cs— Commercial  development 

Industrial  Laud  Area 
Areas: 

It— Tank  farm 

Ig— Generating  station 

Ic— Cement  plant 

Io— Other  manufacturing 

Extractive  Industry 
Areas : 

Es— Stone  quarries 

Eg-Sand,  gravel,  clay,  borrow  pits 

Outdoor  Recreation 
Areas : 

OR-AH  outdoor  recreation  facilities 
Point  data: 

Outdoor  recreation  facilities  (OR):  typespresent 

Golf  courses  (OR-1) 

Ski  areas  and  related  winter  sports  (OR-2) 

Swimming  beaches  and  pools  (OR-3) 

Marinas,  yacht  clubs  (OR-1) 

Stadium  (OR-6) 

Drive-in  theaters  (OR-7) 

Fairgrounds  (OR-8) 

Parks  (OR-9) 

Race  tracks  (OR-11) 

Amusement  parks  (OR-12) 

Zoo,  botanical  gardens,  museums  (OR-M) 

Public  and  Semi-Public  Institutional  Land  Uses 
Areas : 

P— All  areas 
Point  data: 

Public  and  semi-public  areas  (P) :  types  present 
Educational  (P-l) 
Religious  (P-2) 
Health  (P-3) 
Military  (P-4) 
Cemetery  (P-6) 
Water  treatment  (P-7) 
Prison  and  other  correctional  (P-S) 
Governmental  (l'-9) 
Sewage  and  Solid  Waste  Disposal  Facilities 
Areas: 

D-All  areas 
Point  data: 

Sewage  and  solid  waste  disposal  facilities 
(D):  types  present 

Sewage  treatment  plant  (D  I) 
Refuse  incineration  (D-2) 
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Transportation  Laud  Area 

Arc...: 

Th— Highway  interchanges,  limited  access 
right-of-way,  etc. 
"I'r — Rail  yards 
Ta— Airports 
Tb-Barge  Canal  (active) 
Tp— Port,  clock  facilities 
Ts— Shipyards 
Tl-Locks,  dams 

Tt— Communications  and  utility  facilities 
Point  data: 

Highway  category  (h) :  highest  present 
None  (h-0) 
Minor  road  (h-3) 
Two-three  lane  (h-4) 
Four-lane  (h-5) 

Divided  (h-6)  .  .: 

Limited  access  (h-7) 
Hudson  River  bridge  (h-8) 
Major  interchange  (h-9) 
Railway  facilities  (Tr):  types  present 
Abandoned  right-of-way  (r-1) 
Active  track  (r-2) 
Stations  and  structures  (r-4) 
Bridges  (r-6) 
Airport  facilities  (Ta) :  types  present 
Personal  (a-'-t) 
Noncommercial  (a-2) 
Commercial  (a-3) 
Ail  line  (a-4) 
Military  (a-5) 
Heliport  (a-6) 
Seaplane  base  (a-7) 
Utilities  (tt):  type  present 

Electric  power— long-distance  transmission  (t-1) 
Historical  Transportation  Area 
Areas: 

I  lb— Abandoned  Barge  Canal 
Environmental  and  Natural  Resource  Factors 
(SUMMAKY  SHEtT) 
Soils:  predominant  type 
5-1  Organic 
S-2  Clay 
5-3  Silt 

5-4  Sand  and  gravel 
5-5  Till 
S-6  Sand 
Depth  to  bedrock:  predominant  condition 

B-l    o-  ;.'  i 

H-Z  6-10' 
B-3  11-20' 
li-1  21'  t 
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Ground-water  table:  predominant  class 

G-l  High  (surface  to  -10') 

G-2  Low  (below -10') 
Land  forms:  (ypt*s  present 

L-  I  Delta  * 

L-  2  Flood  plain 

I.-  3  Terrace 

I.-  4   Beach  lines 

L-  S  Lakcbcd 

1,-  6  Outwash  plain 

L-  7  Karnes 

L-  8  Eskers 

L-  9  Drumlins 

L-10  Moraine  (ground) 

L-ll  Escarpments  (at  least  "1  km  in  length) 

L-12  Hills  (local  relief  less  than  700') 

L-13  Mountains  (local  relief  700' 4- ) 

L-34  Moraine  (terminal) 
Ridge  Lines:  presence  or  absence 

R-0  not  present  in  cell 

R-l  present  (extending  over  1,5  km) 
Edges:  presence  or  absence 

J-0  not  present  in  cell 

J-l  present  • 
Landmarks:  types  present 

7.-1  Bridges 

Z-2  Chimneys 

2-3  Towers 

Z-4  Steeples 

2-5  Gorges 

2-6  Rock  outcrops 

2-7  Waterfalls 
Building  height:  Dominant  height  in  cell 

H-l  l-2!/3  stones 

H-2  3-6  stories 

H-3   7  +  stories 
Building  coverage:  approximate  percentage  in  cell 

F-l      0-20% 

F-2  2,1-40% 

1-3  41-60% 

F-4  61-80% 

T-5  81-3  00% 
Paved  coverage:  approximate  percentage  in  cell 

C-l     0-20% 

C-2   21-40% 

C-3  41-60% 

C-4  61-80% 

C-5  81-100%        " 
Urban  areas  without  trees:  predominant  condition  in  cell 

U-0  Not  a  predominant  condition  in  cell 

U-l   Predominant  condition  in  cell 
Shoreline  condition 

River  shoreline  with  bulkhead  (k):  miles 

Land,  not  currently  used,  between  railroad  and 

river:  area 

Xa— with  access 
X— without  access 
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From  Supplentfnt.il  Maps 

Elevation:  dominant  elevation  in  cell 

e-  1  0-  500' 

e-  2  500-1000' 
e-  3  1000-1500' 
c-  4  1500-2000' 
c-  5  2000-2500' 
e-  6  2500-3000' 
c-  7  3000-3500' 
e-  8  3500-4000' 
c-  9  4000-4500' 
e-11  4500-5000' 
c  12  5000-5500' 
Bedrock:  predominant  type  by  physiographic  province 
Adirondack  Mountains 

b-  1  Granite-gneiss 

b-  2  Metamorphosed  sediments 

b-  3  Marble 
Hudson  lowlands 

b-  4  Shale 

b-  5  Shale-limestone 

b-  6  Limestone 
Hudson  Highlands 

b-  7  Metamorphosed  igneous 

b-  8  Gneiss 

b-  9  Shale 

b-11   Limestone 
New  England  Upland 

b-12  Granite-gneiss 

b-13  Metamoiphics 

b  14   Marble 
New  Jersey  Lowlands 

b-15  Sedimentary 

b-16  Volcanic  sill 

b-17  Gneiss 
Slope:  area  data  (percentage  of  cell  in  each  class) 
0-  5% 
6-15% 
16-25% 
25%  + 
Historic  Sites: 
number  of  sites  by  types 

Types  are  determined  by  an  alphanumeric  code  con- 
sisting of  a  letter  indicating  significance  and  a  digit 
indicating  date  of  origination  (e.g.,  A6,  C$,  etc.): 
A  National  importance 
13  Slate  importance 
C  Local  importance 
1  1609-1644       Dutch  Settlement 
?.  1665-1776       English  Settlement 

3  1777  17S3       Revolutionary  War 

4  1781-1825       Federal  Period 

5  1826  1916       Industrial  Period 

6  T'l 7  1945       Urbanization 

7  194(3  pi. -cut  Rut  World  War  II 
^  no  dale  indicated 
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